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FOREWORD 



This volume contains concise versions of the papers presented 
at the Conference on Environmental Effects on Antenna Performance 
which did not appear in Vol. I of the Proceedings. The papers pre- 
sented here were reproduced from the manuscripts supplied by the 
authors. Editing on the submitted material v,as relatively light, 
although many of the manuscripts required at least some revision. 
As editor, I should like to thank the members of the advisory board 
for their assistance and, in particular, I am grateful to the con- 
ference secretary, Mrs. Eileen Brackett, for her diligence in 
assembling the papers and readying the material for the printer. 

Volume I which was issued earlier is available from the 
Clearinghouse for Federal Scientific and Technical Information, 
Sills Building, 5285 Port Royal Roa.d, Springfield, Virginia 22151. 
When ordering specify that the accession number is PB-183 715. 
Additional copies of this Vol. II are also available from the CFSTI 
for $3.00 and the accession number is AD 392 879.. 



James R. Wait 
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rROIATE SPHEROIDAL AND LINEAR ANiSKMS IN LOSSY MEDIA 



R. J. Lytle 

Lawrence Radiation Laboratory, University of California 
Livermore, California 



Abstract. The short circuit current and charge distributions and 
the complex effective length for a prolate spheroidal receiving 
antenna, in a complex medium, have been determined for an arbitrary 
angle of incidence. In addition, the power decay as a function of 
§ has been evaluated for a variety of " situations for a prolate 
spheroidal transmitting antenna. 

Tb boundary value problem depicted in Figure 1 has been solved 
by utilizing the method of separation of variables. 
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*Work performed under the auspices of the U. S. Atomic Energy 
Commission. 



<[(§ -1) m \rf cos cos m0 + rag R u; (c,g) sin sin v0j A + 

C °S "in 

? dR ( ^(c,g) a v 
rar) [(s -l) mn 5 g cos cos m0 + mg R^'(c,g) sin sin m0] - 

P dR^(c,5) 7 

a [(g -l) — J5iU sin sin n0 + mg R^'(c,§) cos cos ra?]V, 

mn os mn j 

where 3 and a are the expansion coefficients for the scattered 
field, & fi 4 j» ejW (Cti^jo). 

As E T = at g --• § , ry multiplying (3) by S^(c,7|) and inte- 
grating the results from -1 to +1 in 'fj and from to 2n in 0, and 
evaluating the result at g = g^, one can obtain for E q = k 



p—dR^^Cg ) g , v 
fJe^ i It O , O „(lj/ e M 



o 



By utilizing this result, one can show that the current 
distribution l(1|) for the short circuited receiving antenna is 



(5) i GC (ro --^ UiUr.i-f) ' 



kR 

c 



The short circuit charge distribution Q (T]) can be determined via use 
of the continuity equation and equation (5). 
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By appropriately coupling *.he above results with previously 
obtained results (Lytic and Schultz,1968) for the current distri- 
bution I t,r pi), charge distribution Q tr (T|), and input impedance 2 C 
of a prolate spheroidal transmitting antenna, one can obtain such 
quantities as the complex effective length and the current and 
charge distributions of a loaded receiving antenna. 

Another quantity of interest is the rate of decay in the com- 
plex medium of the power generated by a transmitting antenna. This 
has been determined by evaluating the power passing through a series 
of confocal prolate spheroids surrounding the transmitting antenna. 

Examples of results obtained for some of the above quantities 
follow . 

R. J. Lytl<j and F. V. Schultz (1968), Linear Antennas in Plasma 
Media, Purdue University, School of Electrical Engineering, 
TR-EE 68-10. 




Radiated Power Dependence Upon Phase Length and Loss 

("a" is the seminsinor axis of a confocal spheroid 

surrounding a transmitting antenna with a semimiuor 

axis "a ") 
o 
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Short Circuit Current Dependence Upon Phase Length and Loss 




Length and Loss 
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NUMERICAL STUDIES OF THE EFFECTS OF NON-PLANAR 
LOCAL TERRAIN AND GROUND SCREENS 

by 

R. V. Row and D. M. Cunnold* 
Sylvania Electronic Systems, Waltham, Mass. 

ABSTRACT 

Numerical integration of the vector form of Green's theo- 
rem for electromagnetic fields has been used to evaluate 
the coupling between a pair of vertical dlpoles over a 
two-dimensional non-planar electrically non-homogeneous 
ground. Antennas with arbitrary free-space radiation pat- 
terns can be accommodated. The results of this method are 
compared with previous computation for planar ground and 
ground screen using the compensation theorem approach. 
Agreement is good except at angles below a few degrees. 
New results on the effects of a finite ground screen sys- 
tem in proximity to a sloping beach and a flat sea with 
varying tide levels are given and compared to some HF 
band measurements. 



In an earlier paper [Cunnold, Row, Arnold, 1968] the wri- 
ters described a method of determining approximately the 
coupling between two horizontal electric (dipole) antennas 
over a piecewise homogeneous non-planar two-dimensional 
ground. 

The method if based upon numerical integration of Green's 
vector theorem for electromagnetic fields 

E s - ~ s / [-iwp(nxH)9 + (nxE) x + (n.E)^lds (1) 
e -ikR 

where $ = — r — is the scalar free-space Green's function 



* Dr. Cunnold is a consultant to Sylvania 
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and the required field components nxH, nxE and n«E are ap- 
proximated at each point of the surface in the sense of 
classical optics, but using appropriate Fresnel reflection 
coefficients at each point of the surface. Agreement with 
a limited set of experimental measurements was seen to be 
excellent. Liepa ( 1968] has applied a similar technique 
to determine the field from a horizontal line current 
source illuminating a hemi-cyllndrical boss on an infinite 
conducting plane, a problem for which an exact analytic 
solution is available. He found excellent agreement be- 
tween the exact and approximate field distributions. 

The approach used by Cunnold et al { 1968) has been ex- 
tended to vertically oriented dipole antennas over two- 
dimensional ground, as shown in Figure 1. As in the ear- 
lier work, the ground surface is approximated by a series 
of connected planar strips of infinite extent in the y 
direction. The computer program evaluates the normalized 
transfer impedance 

j^^^K 1 (2) 
rad rad i rad 

where R . is the radiation resistance of an elementary 
dipole, rac zJ[ is the transfer impedance between the anten- 
nas in the absence of the ground, and z\ is the ground 
scattered contribution from the i-th strip. 

For a short electric dipole of length df 

Z = E dl 
t x 

Equation (1) is used to evaluate E x so that after a sta- 
tionary phase integration in the y f coordinate only the z' 
integration need be done numerically. 

inc 

The field components, H«t, E^,, E,i and E x are calcu- 
lated on the basis of Norton' s [ 1937) far-field expres- 
sions. 
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Comparlson with Known Results 

Provided the Norton ground wave terms are included in the 
field expressions and the lower antenna of the pair is 
taken as the receiver, excellent agreement is obtained be- 
tween the strip method and Norton's results for dlpoles 
over a homogeneous flat ground. Wait [1967} has examined 
the situation of an elevated electric dipole over a cir- 
cular ground screen on a flat earth. Agreement with Walt's 
results is excellent except at angles below 5°, where the 
non-stationary nature of Eq. (1) will exhibit errors in 
the final result through errors in approximating the true 
fields nxH, nxE and n«E to a greater degree than the sta- 
tionary formulation used by Wait. 

Effects of Proximity to the Sea 

An antenna installed close to the sea is in an environment 
which changes with tide level. Figure 2 shows a profile 
for which extensive computations have been made at fre- 
quencies of 5, 10, 2,0 and 35 MHz and screen lengths of 
150, 500 and 750 feet. The receiving antenna is 1 foot 
above the screen, representing a monopole. Figure 3 shows 
the effect of increasing the screen length from 150 feet 
to 500 feet to 750 feet under low tide conditions. A sig- 
nificant improvement in gain is realized in going from a 
150- foot to 500-foot screen, with the maxixsum improvement 
occurring at a progressively lower angle as the frequency 
is increased. A lesser additional improvement in gain is 
realized by extending the screen over the down sloping 
section of beach to the edge of the sea. Figure 4 shows 
the effects of a 10- foot increase in tide level with a 
150- foot screen. The maximum change (increase) in gain 
is about 1 dB and shows up at the highest frequency. Ad- 
ditional computational results not presented here show 
that changes in gain caused by tide level are somewhat 
smaller for the greater ground screen lengths. 

The computations shown in Figures 2 and 4 were repeated 
for a receiving antenna with a cos 6 pattern. Apart 
from a general reduction in power gain at angles above 
about 10°, the results are in substantial agreement with 
those shown in these which, of course, pertain to a dipole. 
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Comparison with Observation of Tide Level Changes 

Extensive flight test pattern measurements [ Simons et al 
1968] have been made on an HF band log-periodic monopole 
element installed on a low profile beach installation at 
Keflavik Naval Station, Iceland. The contour lines are 
more or less parallel over the 260° to 340° .seaward sec- 
tor and hence constitute a reasonable approximation to a 
two-dimensional ground over this sector. The most in- 
teresting of these data for our purposes (shown in Figure 
7) are the 30-MHz measurements taken on October 29, 1967, 
at an azimuth of 260° (which extends over the sea). One 
set of these dcta pertain to low tide and the other to 
high tide (water 8-feet higher). Figure 5 shows the 
actual profile ( solid curve) along the 260° azimuth, 
where the datum le\el is transferred from mean low water 
(low tide) to tne 14 -foot contour for purposes of enter- 
ing the baseline elevation data into the computer pro- 
gram. The lightly dashed straight line is the sloping 
beach profile approximation which simplifies the labor 
of entering profile data into the computer. The strip 
method program was run for both these profiles; the 
former (actual profile) case being designated strip-by- 
strip format in Figure 6. In the computations the 
ground screen (No. 10 AWG wire with 12-inch spacing) ex- 
tended from -55 feet to +80 feet. It is evident from 
Figure 6 that the results obtained with the sloping 
beach format are in excellent agreement with the strip- 
by-strip input, especially for the low tide case. The 
additional set of curves, shown with light dash and 
double dots in Figure 6, refer to the gain of a pair of 
vertical dipoles 6 feet and 14 feet above an infinite 
flat sea, respectively. These latter curves show that 
in general it is not possible to predict the influence 
of changes in sea level by such a simple model. There 
is a striking resemblance between the detailed shape of 
the predicted curves as shown in Figure 6 and the mea- 
surements including a crossover at about 14° in both 
measurements and predictions. The maximum increase in 
gain in going from low tide to high tide is 2.1 dB at 
an angle of 10° in the predictions, and 3.6 dB at 8.5° 
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in the measurements. It should be noted that at angles 
below 1.5° the high tide measurements show a relative 
field strength that is about 1.4 dB greater than at low 
tide. The opinion is offered that these two relative 
field strength (gain) curves should be coincident at all 
angles below 1.5°, as is the case for dipoles at low 
heights over an infinite sea shown in the dash-double* 
dot curves of Figure 6. The dash-double-dot curve in 
Figure 1 represents a uniform downward shift of 1.4 dB 
for the high tide experimental data. When this is done 
the maximum "observed" increase in gain in going from 
low tide to high tide is 2.4 dB at 9°, which compares 
even more favorably with the strip method predictions in 
Figure 6. Of course the "observed" crossover point is 
now shifted down to 12.5° compared to about 14° in the 
predictions. 

It is concluded that the strip method does very well in 
predicting the detailed shape of the radiation pattern 
at angles above 1.5° or 2° in this one instance of a 
low-profile beach installation. The fact that the mea- 
surements were made using the LP element on the beach 
whereas the predictions are for a short vertical dipole 
should not change the conclusions substantially, as evi- 
denced by the discussion in the preceding section. 
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Figure 1. Coordinate Systems for Computing Transfer 
Impedance 
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Figure 2. Geometry of Merrimac Light Coast Guard Station Sloping 
Beach Profile 
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THE CURRENT DISTRIBUTION ON A FINITE LENGTH 
DIPOLE IN THE PRESENCE OF GROUND SCREENS 

Virgil R. Arens, U. R. Embry, and D. L< Motte 

Sylvania Electronic Systems, Western Division 
Post Office Bet 205 
Mountain View, California 94040 

ABSTRACT 

Some of the mathematical and numerical techniques 
using digital computers which will lead to determining the current 
distribution on a finite length dipole in the presence of radial 
ground screen are presented. Computed results for the input 
impedance for non-resonant vertical dipoles near tne interface 
for various ground conductivities and dielectric constants are 
given. The input impedance for a non-resonant dipole in the pres - 
ence of a planar array of rad.lal wires of various lengths is shown. 

1. INTRODUCTION 

On a number of occasions investigators {Wait, 1967) 
have had to make assumptions with respect to the current dis- 
tribution on finite length dipoles above and near variously shaped 
ground screens. It is our goal to develop a set of techniques 
that will allow one to calculate, using digital computers, the 
current distribution under these conditions. 

The methods under investigation are the method of 
moments (Harrington, 1968), Hallen's integral equation (Hallen, 
1938; Mei, 1965), and the integral equation approach (Tanner 
and Andreas en, 1967). All three of these approaches make use 
of the integral of free space Green's function. 

. 2. Green's Function 

In our study of dipoles, we believe that we have devel- 
oped some useful methods in evaluating the integral of the Green's 
function. 

In the first method one can simply perform a numerical 
integration of the function. Perhaps one should not use the word 
"simply", though, because the integration for the diagonal terms 
of the matrix pose some extremely difficult numerical problems. 
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For these terms the observation point lies on the surface of the 
integration segment so that the integrand comes within $ of a 
pole, where 6 i.s the radius of the wire. To get an idea of the 
difficulty this introduces, the approximated numerical value for 
the integral has not become stable after using up to 64 points of 
Gaussian quadrature (corresponding to a polynomial of degree 127 ) 
and up to 9000 points -.sing the trapezoid rule. 

In the second method one rewrites the equations in cylin- 
drical coordinates. This considerably simplifies the work because 
p and z components can be integrated separately after a local rota- 
tion and translation to give a two dimensional problem and then 
dotted with a direction vector at the observation segment. 

The third method makes iwe of spline functions an* 1 has 
the potential of allowing us to develop a program free of all nu- 
merical integrations. This is accomplished by a piece -wise 
approximation of e -jkR in cubic polynomials over one cycle. When 
the resulting polynomials in the va riable kR are placed under the 
integral it is found that the terms ox the resulting integrand are all 
integrable in closed form. The particular functions one uses are 
called cubic splines, as excellent account of which is found in 
Hands comb (1966). 

As indicated earlier, when the observation point lies on 
the segment of the cylinder that is to be integrated, a singularity 
occurs in the expression for the electric field. In point of fact in 
the derivation of the equation an interchange of differentiation and 
integration was made which is not necessarily valid for these seg- 
ments. Fortunately, VanBladel (1961) investigated this case some 
years ago and h?d defined a principal value of the integral for the 
diagonal terms. 

3. Evaluation of the Integral of the Grqen's Function 



The evaluation of the integral of the Green's function as 
was indicated earlier can at times be evaluated by using only 
numerical techniques, or a MacLaurin series expansion. Some 
authors have in the past tried more difficult methods that involve, 
for example, Hankel or Besscl functions. Now note that: 




a a. 
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With this simple change, the integral in the final equation can now 
be accurately evaluated by as few as two or three points of Gaus - 
sian quadrature for all values of R and accuracy of 5 significant 
digits , since the integrand is finite for R equal to zero. Here, 
then, is a simple and accurate numerical technique for evaluating 
the f reespace Green's function, which we have not seen elsewhere 
in the literature. It illustrates very nicely the fact that it is not 
always expedient to begin a numerical evaluation of a function by 
expanding it in a power series. 

4. Antennas Over A Finite Conducting Ground 

■ i i i iumiiiwi — — ■ i i i i mi mm - l inn i n m i i Hh— iii w mimii — 

In considering antennas over ground, it is easy to see 
that the free -space Green's function is no longer valid. Here one 
needs to match a set of boundary conditions at the interface be - 
tween the half spaces. A Green's function which must meet these 
conditions is called a modified Green's function. Banos (1966) 
gives the properly modified Green's functions that are valid nearly 
everywhere in the upper half space. 

The total Green's function is the sum of the free space 
Green's function for the image of the source segment and a 
correction term to account for the ii.zt that the ground has finite 
conductivity. Banos deals primarily with the evaluation of this 
correction term. If one considers a dipole just above the ground, 
then one finds that none of the expressions given by Banos are 
valid. For this region we are developing an adaptive numerical 
integration technique. 

5. Ground Screens 

Ground screens play two roles in high frequency antennas. 
The first role is used in association with monopoles and the 
ground screen then is commonly called a cc iter poise. Here one 
is interested in distributing current into the earth in a least loss 
manner to give an antenna with a high efficiency. The second 
role, which is used with both monopoles and dipoles, is the use 
of ground screens to modify or change the far field pattern of 
antennas. 

The method of obtaining the current distribution in the 
presence of ground screens can be obtained in at least two ways. 
THe desired method is dependent on the number of wires in the 
ground screen. Consider for a moment, the simplest of all counter- 
poises, a limited number of radial wires buried in the earth used in 
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conjunction with a monopole. In this system the counterpoise is 
considered a part of the antenna. The impedance matrix which 
defines the antenna includes terms representing the mutual im- 
pedance between segments of the antenna and the buried wires. 

When the number of wires in the ground is large, the 
size of the matrix may become too large to consider the buried 
wires as part of the radioing system. When the size of the 
ground screen is sufficiently large one can take the expressions 
for surface impedance given by Wait {1956, 1958) and Lars en 
(1962), and determine an equivalent dielectric constant and earth 
conductivity and substitute these in Bano's equations to determine 
the modified Green's function. 

If the ground screen is not large or the edge of the screen 
is near the antenna it is conceivable that one could determine a 
different modified Green's function. However, in view of the fact 
that one is attempting to determine the mutual impedance between 
two Hertzian dipoles due to the presence of ground screens. 

6. Reducing the Number of Equations 

Use of the symmetries of a system reduces the work in- 
volved in solving problems connected with the system. In simple 
structures, the symmetry conditions may be clear, and intuitive 
ideas are adequate to take advantage of them. In a more compli- 
cated system, intuition may become vague and a systematic method 
is needed to examine the effects of the various symmetry condi- 
tions. We have developed a two step procedure that considers all 
b/mmetxy conditions when setting up the calculation of the current 
at specified points where the current will be determined. Symme- 
try conditions may show that some of these points will have equal 
currents; hence, assign only one variable and gather the coeffi- 
cients to form a single term in each equation. (This is equivalent 
to adding columns of the equation matrix) This leads to an over- 
deter mined system, with more equations than unknowns. The 
second step is to determine which equations are redundant and 
remove them; thus obtaining smaller sets of equations to be solved. 

As an example, consider a monopole with radial ground 
wires spaced ten degrees apart. If one specifies that the current 
will be found at five points on each wire, the result is a set of 
185 equations in the 185 unknown currents. If the symmetry of the 
ground screen is invoked, the current distribution in any wire must 
be the same as any other, so five unknowns suffice for the ground 
screen current values. The set now involves ten equations in ten 
unknowns. 



-25- 



7. Results 

The current distributions and input impedance were com- 
puted for a center -fed -non -resonant dipole. The dipole length 
was 15 meters, the diameter was 0. 006 meters and the fre- 
quency was 10 mHz, For compactness, only imput impedance plots 
are shown. For most cases, the current distribution obtained was 
not sinusoidal. 

Figure 1 shows the input impedance oi the dipole in free 
space, and with the tip 0, 01 meter above the surface, in a vertical 
position, over various types of ground. The value fpr the perfect 
conductor was checked by computing the impedance for each of two 
collinear antennas fed in phase. 

Figure 2 shows the effect on the input impedance when a 
radial array of wires is placed near one end of the dipole. Values 
are shown for various lengths of 4 -wire and 3 6 -wire arrays. The 
values for 9 -wire and 18 -wire arrays show similar patterns. As 
might be expected, many short wires have more effect than a few 
short wires. As the wires in the array approach a length of half 
a wave-length, the arrays with few wires cause a greater change 
than the arrays with many wires. 

Placing a vertical dipole over smooth homogenous ground 
appears to increase the input resistance. Placing the dipole over 
an array of short radial wires increases the reactance, with a 
small change in the resistance. The combined effects are ex- 
pected to cause the input impedance to approach that of the dipole 
over a perfect reflecting surface. K the array contains many 
wires, the length should be about 0. 1 5 to 0. 2 wave-length; with 
few wires, the length should be greater, to about 0. 25 wave length. 

9. Conclusion 

The digital computer programs that are being developed 
will, when complete, allow one to evaluate the properties of any 
dipole antenna in any orientation between horizontal and vertical 
over a ground screen. 

The approaches are being developed with the idea in mind 
that other antennas such as the log periodic or rhombic are more 
important than dipoles ; the work on dipoles is a means to develop 
the techniques. 

Besides the advantage of speed, economy, and flexibility 
relative to experimental methods, the methods outlined here will 
provide answers with a precision and detail unobtainable by flight 
testing full scale high frequency antennas. The knowledge of the 
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current distribution and input impedance provides insight to the 
antenna's operation and enhances the engineer's design tools. 

REFERENCES 



Wait, J.R. (1967), Canadian Journal of Physics. Vol 45, p. 3091, Sep 1967. 

Harrington, R. F. (1968), Field Computation by Moment Methods. 

MacMUlian Co, New York, 1968. 
Hallen, E. (1938), "Theoretical Investigation into the Transmitting and 

Receiving Qualities of Antennas, " Nova Acta. Royal Soc. Sciences . 

( Uppsala) . Vol 11, pp 1-44, Nov 1938>x 
Mei, K. K. (1965), "Onthe Integral Equation of Thin Wire Antennas," 

IEEE Transactions on Antennas and Propagation . Vol AP-13, 

pp 374-378 (May 1965). 
Tanner, R. and Andreas en, M. (1967), "Numerical Solution of Electro- 
magnetic Problems, " ffiEJE_Sp_ectrum, Vol 4, No. 9, pp 53-61, 

Sept. 1967. 

Silver, S. (1949). Microwave Antenna Theory and Design . McGraw-Hill 

Book Company, 1949. 
Van Bladel, J. (1961), "Some Remarks on Green's Dyadic for Infinite 

Space, IEEE Transactions on Antennas and Propagation, Vol AP-9, 

Nov 1961, pp 563-566. 
Handscomb, D. C. (1966), Methods of Numerical Approximation . Pergamon 

Press, pl63, 1966. 
Banos, A. (1966), "Dipole Radiation in the Presence of a Conducting Half- 
Space, Pergamon Press, 1966. 
Wait, J.R. (1956), "On the Theroy of Reflection from a Wire Grid Parallel 

to an Interface Between Homogeneous Media I, " Applied Science 

Research , Vol B6, p259, 1956. 
Wait, J. R. (1958), "On the Theory of Reflection from a Wire Grid Parallel 

to an Interface Between Homogeneous Media H, " Applied Science 

Research . Vol B7, p 355, 1958. 
JLarsen, T. (1962), "Numerical Investigations of the Equivalent Impedance 

of a Wire Grid Parrallel to the Interface Between Two' Media, " NBS 

Jc urnal of Research, Section D, Radio Propagation, Vol. 66 D, No. 1 

pi\ 1962. 

Wait, J. R, "The Theory of an Antenna Over an Inhomogeneous Ground 

Plane, " Electromagnetic Theory and Antennas, (Ed. E. C. Jordan), 
Pergamon Press, Oxford, pp. 1079-1098; 1963. A 

Wait, J. R. and L. C. Walters, "Influence of a Sector Grounc" Screen on the 
Field of a Vertical Antenna* NBS Monograph 60, April 1963. B 



-27- 



§ 

W 

u 

h 

< 

W 

en 



80 



60 





— i r 


1 r~- ■■« 


1 




J" 










e 


a mhos /m 








1. 


1.0 


0.0 








2. 


3.0 


0. 001 








3. 


5.0 


0.005 


- 






4. 


10.0 


0.02 








5. 


80.0 


5.0 








6. 


1.0 


1.0 x 10 7 


■ 


1 












— i — ...i . 


—i 1 X. 




i 


L ., 



to 



tOO 120 140 

RESISTANCE, OHMS 



Fig. 1. Input Impedance of a Half- Wave Dipole Over Ground 
lOOT 




80 100 ' "120 MO l — 

~. , T 4 T „ , RESISTANCE, OHMS 

Fig. 2. Input Impedance of a Half- Wave Dipole Over a Radial 
Array of Wires of Varying Length. 



-28- 



MEASURED PATTERNS OF HF ANTENNAS AND 
CORRELATION WITH SURROUNDING TERRAIN 

D. R. McCoy, R. D. Wengenroth 
General Electric Company, HMES, Syracuse, N. Y. 

and 

J. J. Simons, RADC, Griff iss Air Force Base, Rome, N. Y. 

A large number of fullscale radiation patterns have been measured on several HF 
antennas with various siting conditions. From the data it can be demonstrated conclusively 
that siting plays a major role in the radiation characteristics, as expected. A study of the 
data reveals ways of modeling the earth so as to adeqvattly predict the influence of various 
parameters on antenna coverage. Thus, these data provide a firm reference with which 
various theories can be verified. Only two examples from the many available are presented 
in this paper. 

1. Utroductioa 

It has been well established by theory that antenna surroundings will modify the far-field radiation patterns 
of an antenna system -that in fact, the antenna is composed of the actual radiators and any territory many miles 
away which may be illuminated by the system. In the high frequency range, where the wavelengths are quite 
large and the antenna beam widths in general relatively broad, this phenomena becomes extensive. This paper 
presents some results of a measurement program on such a system, demonstrating the effect Only a brief 
sampling of the data is presented here. A more complete presentation is made elsewhere/ 2 '^) It also discusses 
techniques whereby predictions may be made of the impact of surroundings on such a proposed installation. 

2. Discussion 

The full-scale patterns were measured using a KC-135 aircraft, instrumented with a receiving antenna on the 
tail boom, and a set of four receivers multiplexed inside. Four transmitters operating at different frequencies 
were time-multiplexed into the antenna on the ground, allowing for four patterns to be measured on each pass 
over the antenna. The aircraft was flown a! a constant altitude (nominally 30,000 ft.) and the range and signal 
strength were recorded in the air. Flight path was controlled by ground radar and the navigator out to a 
maximum tange of approximately 2S0 nautical miles. The data was then scaled, calibration data applied and the 
elevation angle and signal strength calculated and plotted, yielding relative data which in most cases was 
reproducible to within I dB. 

These full-scale patterns of HF antennas installed with various siting conditions have been made jointly by 
General Electric Company and the RADC test center The installations range from a site atop a mountain 
overlooking the ocean in Greenland to a sea shore site in Iceland, with the ground screen end touching the ocean 
water at high tide. Most of the results discussed here are of two identical sets of log periodic elements, one 
looking north and the other south over land in upper New York State. Of particular interest is the series of 
patterns obtained at boresight (188°) and at plus and minus 10° (198°, 178°) of boresight looking south. The 
antenna is on a hill overlooking the Mohawk valley. From the elevation of 1760 feet above sea level the land 
drops off rapidly to 7-900 feet (a slope of approximately 6°) into the valley 10 to 20 miles away. In the 198° 
direction the terrain is broken up with several hills so that there is no flat land. This pattern (not shown) does 
not display the rapid lobing at the higher angles. The profile shown in figure 1 demonstrates the condition on 
boresight (188°). The radiation pattern for this direction is shown in figure 2. The identical model of antenna 
looking to the north from the same hill exhibits entirely different radiation characteristics, as demonstrated in 
figure 3. In this direction the initial slope is about 6° but the elevation levels off at about 1400 to 1500 feet in a 
short distance and stays constant to the Adirondack Mountain region 30-40 miles away. This contouring is 
shown in figure 4. 

A consideration of the specular reflection from a combination of the slope in close and the relatively flat 
terrain further out leads to a rapidly varying pattern at low angles and a wide minimum at sli°h"y higher 
elevation angles. It therefore appears from an observation of measured data that there is a third effect at the 
lower angles, in which the flat plane ripple is suppressed. This is probably due to the energy being blocked by the 
slight roll immediately in front of the antenna. There were minor differences in the lobing structure, particularly 
to the south, when snow was on the ground. Several values of ground constants were considered in calculating 
the above mentioned patterns. Proper choice of ground constants (K = 30 and 0= .001) seemed to give a good 
fit to the measured data, with the reservations of suppression of rapid lobing at the lower angles. 
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Figurs 5 shcvs the measured patterns ai ihc Kcfiavik (iceiand) site. The anienna is located on the shore line 
with the ground plane extending to the water. In this case there was a noticeable effect of the tidal variations, 
which were greater than 13 feet. A calculation of the expected patterns for the low and high tide conditions 
provides an exccUent prediction of the effect as can be seen from figure 6. This same antenna, when set on a 
2S60 foot hill overlooking the snow, rock, hills etc. and finally the occwi 17 miles away, at Thule, Greenland, 
displayed violent lobing, even though the earth dropped off sharply in fret of the antenna to an elevation of 
approximately 1500 feet above sea level. 




-•.0 -40 4.0 CO ••0 -4.0 4.0 CO 
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Figure 5 Figure 6 

A technique of modeling the land in front of an antenna with a scries of flat planes and knife edges has been 
developed/ This program estimates the worst case limit for the particular geometry assumed. A result of using 
this program to analyze the Thule antenna is shown in figure 7. In this case the very rough terrain was assumed 
to be approximated by two knife-edge points and intervening flat planes. As can be seen, the low angle coverage 
prediction is quite good. It is obvious that an added effect occurs in the 2.5° to 4.5° region, causing a wide lobe 
of high amplitude and phasing with the more regular pattern. Perhaps a superposition of these data with returns 
from additional specular points would be in order to further refine this result. 

In the case of the data presented above for the Starr Hill site in New York State, there are no significant 
sharp obstructions in the path to the north and the program does not apply; however, it has been used with fair 
success on the south-looking antenna/ 2 ) The program utility can be enhanced by more accurate modeling and 
perhaps by replacing the knife edge diffraction with hemicyiinders and allowing for extended specular reflection 
planes. This is being studied at the present time. 
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RELATIVE GAIN (dB) 
Figure 7 



3. Conclusion* 



From the observed data it is completely demonstrated that siting plays a major part in determining the 
radiation characteristics of an antenna working in at the HF frequency range. It is also evident that, by carefully 
studying the terrain condition and contour for several miles in front of the antenna, a good prediction of the 
performance can be made. 
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RADIATION BY A VHF DIPQLE-TYPE ANTENNA 
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ABSTRACT 



This paper describes the performance of a VHF dipole antenna 
located at the tip of an otherwise pointed reentry nosecone. Trans- 
mission and reflection measurements made during a recent flight ex- 
periment are presented and compared with calculations made for an 
appropriate theoretical model. 



INTRODUCTION 

The performance of a thin (k a = .092) VHF dipole-type antenna 
located at the tip of a 9© nose cone during an almost vertical re- 
entry at a relatively constant velocity of 16.8 kft/sec is described. 
The reentry period reported covers a time span of eight seconds with 
an altitude range of 230 to 100 kft. 

The antenna as shown in Fig. 1 is an asymmetrical dipole with 
the cylindrical quarter-wave stub acting as one arm and the payload 
skin as the other arm. It is balun-fed, the outer conductor of the 
cable being connected to the cone and the center conductor to the 
cylinder. A quarter-wave choke section was introduced behind the 
conical arm, hopefully to decouple the aft portion of the conical 
payload. That this was not fully achieved is evident from the 
free-space radiation pattern of Fig. 2. Note the "notch" at broad- 
side which is primarily attributable to the excitation of the 
afterbody. 



*The work reported in this abstract was supported in part by Contract 
F33615-67-C-1643 between Air Force Avionics Laboratory, Air Force 
Systems Command, Wright-Patterson Air Force Base, Ohio and The Ohio 
State University Research Foundation. 
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THE FLIGHT DATA 

The altitude, velocity and angular coordinates of the trajectory 
were obtained by four tracking radars, 2750 to 5500 MHz. The X- 
band telemetry signals were recorded by three receiving stations. 
The VHF transmission was received by four slightly separated CP 
antennas and provided separate VHF-AGC records. 

The telemetered signals of an on-board gyro and accelerometer 
together with known initial reentry conditions and the precession 
frequency provided angle-of -attack data so that the aspect angle 
(the angle between the cone axis and receiver direction) could be 
computed. The ,jreces?ion rate was verified by the cyclic variation 
of the VHF transmission level due to the "notch" in the pattern. 

During the flight the impedance of the antenna was measured 
with a four-probe reflectometer and telemetry recorded as shown 
in Fig. 3. Except for one short time interval, preliminary to 
maximum signal attenuation, the modulus of the reflection coeffi- 
cient was relatively constant. In addition the antenna was operated 
in a receiving mode, a radiometer being used to record the effective 
noise temperature of the plasma at the antenna terminals. This data 
is presented in Fig, 4 where the "average" noise temperature is 
approximately 3500OK. Outside of the time period shown the antenna 
receiver was saturated by VHF interference, so at the time the VHF 
transmission was undergoing maximum attenuation the plasma noise and 
•interfering signal contribution could not be resolved. 

In Fig. 5 the VHF transmitted signal which was recovered from 
the calibrated AGC records of three different receivers is pre- 
sented. The solid curve connects the average received values. A 
maximum signal attenuation of 17 dB is observed at 193 kft. 



THE ANALYTICAL MODEL OF THE ANTENNA 



We approximate the flight situation with a gap-excited finite 
dipole in an infinite, uniform plasma column. The far-zone result 
is similar to that for the Hertzian dipole in the infinite column, 
with the exception of the current-plasma interaction. Using the 
result of Tingl et al . The ratio of the power with and without 
plasma radiated in the direction is given by 
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Fig. 3. ANTENNA 'MPE0ANCE 
DURING REENTRY. 



Pig. 4. EFFECTIVE PLASMA NOISE TEMPERATURE 
AT ANTENNA TERMINALS. 
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where 

(2 ) e . i . _J> 

p 1 + i v/u 

21 is the dipole length (approximately x/2), a is the antenna radius 
and b is the radius of the plasma column. The other symbols and 
functions are the usual. 



THE INVISCID FLOW FIELDS 

Some estimate of the parameter which characterizes the cold 
plasma must now be found. The simplified stream-tube method of 
Seiff and Whiting, 2 which assumes 

(1) isentropic flow along streamlines, (2) flow parallel to 
vehicle axis and (3) the shape of the pressure distribution normal 
to the body is predictable (by first-order blast wave theory). In 
our trajectory, conditions (1) and (3) are probably valid but (2) is 
never satisfied. The angle of attack is small, so we hesitatingly 
ignore the crossflow. We also assume an equilibrium flow for all 
altitudes, which is not exactly the case. 

Using the profiles calculated from the inviscid flow, we find 
a monotonically increasing attenuation from 230 to 100 kft. (A 
traveling wave form of current is assumed on the antenna.) Referring 
to the flight data of Fig. 5 we see this is not the case. 



THE VISCOUS BOUNDARY LAYER 

Maximum attenuation at 193 kft. indicates a shock-boundary 
layer interaction. Due to the question as to the nature of the flow 
and the great boundary layer thickness at this altitude, a quanti- 
tative calculation is subject to question. Nevertheless, we proceed 
with the phenomenological approach to the laminar thermal boundary 
as given by Schlichting3 since some estimate of its effect is pref- 
erable to none. 

The quantity that controls the thermal boundary layer overshoot 
is the Eckert number. We assume a constant temperature (300OK) cyl- 
inder. The inviscid flow calculations indicate that the kinetic 
temperature at the surface remains approximately constant at 3000°K 
throughout the trajectory. Since the velocity is constant then the 
Eckert number is primarily dependent on the specific heat v which 
for air is a rapidly varying function of density in the 2500-5000°K 
range. 4 To effect our calculations we assume that the specific heat 
of the boundary layer is specified by the inviscid flow body stream- 
line. (A rather severe approximation due to the variation of c p 
through the boundary lay3r.) 
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Using the matching point scheme of Evans and Huber, the joined 
boundary layer - inviscid solution at 140 kft., shown in Fig. 6, is 
representative of the results obtained at other stations along the 
cylinder and altitudes. Also shown in rig. 6 is the tu D profile ob- 
tained from equilibrium conditions. 



DISCUSSION 

From other work based on plane wave theory it was found that 
thin symmetrical 1 nhomoreneous layers could be represented quite 
well by equivalent uniform layers. It was observed that the 
equivalent thickness was well approximated by the distance between 
the two points where u) p % .7 w^^. 

With the v and a- found from the boundary layei- results and the 
equivalent plasma radius determined as just indicated, Eqs. (1) and 
(2) were evaluated for eight different altitude points. The atten- 
uation due to the boundary layer is shewn in Fig. 5. To effect better 
agreement with experiment, one must add the attenuation due to the 
shock layer ranging from 0.2 dB at 200 kft. to 4.1 dB at 120 kft. A 
more exact result would be given by a modification of Eq. (1) to in- 
clude two coaxial layers. 

In conclusion, it was found that the VHF transmission by a 
dipole-type antenna imbedded in an electrically thin axisymmetric 
plasma sheath is explained by the viscous boundary layer. Although 
the flow is quite complicated, a constant temperature cylinder is 
assumed, and the dipole is unbalanced, it is believed that the flight 
data and analysis tend to show that the antenna currents are somewhat 
isolated from the plasma by the cool leading edge of the boundary layer 
and the presence of a dielectric ablative coating over a significant 
portion of the antenna. This conclusion is consistent with three experi- 
mental observations: 

1. Only as the greatest attenuation occurred did the antenna 
resonate . 

2. There was no significant pattern change as evidenced by the 
ubiquitous notch. 

3. The driving point impedance (Fig. 3) does not indicate 
tight coupling. 
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EFFECTS OF ELECTRON ACOUSTIC WAVES ON A CYLINDRICAL 
DIPOLE TYPE RF MAGNETO- PLASMA PROBE 

By 

Hiroshi Oya* 
Laboratory for Space Sciences 
Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

The effect of electron acoustic wave generation has been in- 
vestigated for a short cylindrical dipole antenna in an aniso- 
tropic (magneto- active) warm plasma, considering the collisions 
between electrons and heavy particles. Calculations have been 
carried out for a homogeneous plasma using the hydrodynamic ap- 
proximation and the quasi static approximation. It is obtained 
that the resistance due to electron acoustic wave increases 
steeply near the upper hybrid resonance frequency. The theory 
developed in this paper can explain the experimental results 
around the above frequency range in terms of emission of electron 
acoustic waves and electron collisions. 



INTRODUCTORY REMARKS 

The impedance of an antenna immersed in a warm plasma has 
been investigated extensively in the recent decade. Considerable 
progress has been made in the case of the warm isotropic plasma 
(a bibliography will be cited in a full version of this work) . In 
the case of the magneto-active plasma, however, the theory to 
explain the electron acoustic wave effects on RF plasma probe 
experiments (Oya and Obayashi, 1966, 1967; Stone et al 1966), 
which has been carried out in the ionospheric plasma, is not 
adequately developed. 

THEORETICAL PROCEDURE 

In this paper the problem is investigated from the standpoint 
of an RF probe. The short cylindrical dipole antenna impedance 
is calculated in a magnet oactive warm plasma; the EM field is 
neglected since the problem is restricted in a region where the 
electrostatic field is completely dominant. The existence of the 
ion sheath is also neglected since the frequency range of interest 
is near the electron plasma and upper hybrid resonance frequency. 
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The macroscopic hydro-dynamic equation is used to descrive 
electron motions considering the electron collisions between 
heavy particles. Maxwell's equation, the equation of continuity 
and the equation of the state are used being transformed into 
Fourier space with respect to the three dimensional problems. 
All time functions are expressed by the factor eJot where u> and 
t are the angular frequency and time respectively. In this case, 
whole set of basic equations is expressed in a comprehensive form 
by introducing a dielectric matrix which has elements expressed 
as functions of co, wave vector 5, plasma frequency r p . gyro- 
frequency fg and temperature T. 

The basic equation is solved to obtain the potential function 
i as the first step of impedance calculations; in this case quasi 
static approximation is employed. One of important boundary 
condition to solve the differential equation is charge q s , and 
current J s distributions on the antenna. It is assumed that 
the particles in the plasma are completely reflected back through 
a contact region between the plasma and the RF probe, which, in 
Fourier transformed space, effectively gives -jk.J s + jwjijL = 
on the antenna. Hence, once we find the charge function q s , the 
problem can be solved by obtaining the inverse Fourier transfor- 
mation. In this work, q s is assumed to be distributed uniformly 
on a short cylindrical dipole whose axis is parallel to the DC 
magnetic field. This is equivalent to the model of Balmain 
[1965], which is used to obtain the impedance in the case of an 
isotropic plasma; and the same technique is nlso used to obtain 
the antenna impedance by using the e.m.f . method. 

The impedance value is finally expressed by au integration 
in which an integrant is a complicated complex function. This 
integration is carried out numerically for the ionospheric plasma 
parameters . 

RESULTS AND DISCUSSIONS 

A numerical calculation of the impedance value is obtained 
for a short cylindrical dipole antenna whose elements have length 
L -= 1m, and radius pj_ = 0.01m. In this result, the gyroresonance , 
the plasma resonance, and the upper hybrid resonance are indicated. 
The remarkable evidence of the electron thermal effect is the 
presence of the acoustic wave resistance which is enhanced near 
the upper hybrid resonance frequency and, also, in a frequency 
range below t he gyrores onar ■» frequency. In a frequency range 

higher than J2 f 2 + f 2 the electron acoustic wave cannot exist 

p H ' 

since the wave length becomes smaller than the Debye shielding 
length; thus the acoustic resistance is not present. 
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In practice, we cannot neglect the existence of an ion 
sheath surrounding the RF probe, except for the case of sheath 
collaps condition which can be artificially produced by applying 
a high DC voltage to the probe. In the frequency range in which 
|Z(«J | » IVCgiD i where C s is the so called sheath impedance, 
however, the effect of sheath impedance is completely negligible 
[Oya, 1967]. Thus we can confirm that above the sheath resonance 
frequency the ion sheath effect is not serious, and around the 
plasma, and upper hybrid resonance frequency the ion sheath effect 
is completely negligible. 

A microscopic approach to electron acoustic waves in an 
anisotropic warm plasma was considered by Bernstein [1958]. 
Bernstein's relation is developed for a case near the upper hybrid 
resonance by Fejer and Calvert [1964]. The dispersion equation 
is only slightly different from Fejer and Calvert's equation near 
the upper hybrid resonance frequency and indicate that the theo- 
retical approach of this work in a frequency range around the 
upper hybrid resonance obtains results similar to the microscopic 
treatment of the plasma and is sufficient to discuss the quantita- 
tive nature of the upper hybrid resonance. 

In Fig. 1, examples of impedance value using absolute value 
and phase angle are illustrated for plasma parameters of the 
ionosphere at altitude of 100 km (a) , 500 km(b) and 1000 km(c) . 
The phase variation indicates that near the upper hybrid resonance 
there is strong enhancement of the electron acoustic wave; the 
power dissipation of this acoustic waves affects the sharpness 
of the resonance. In the collision dominant regime (see a), 
however, the power dissipation is caused by the collisions. The 
sharpness of the resonance is investigated for the ionospheric 
plasma in altitude range of 100 km ~ 1000 km being compared with 
rocket observation result [Qya, 1967]. There is good agreement 
between the theoretical calculation and the experimental result. 

CONCLUDING REMARKS 

The theory is developed in an anisotropic warm plasma for a 
short cylindrical dipole RF probe immersed in an homogeneous 
ionospheric plasma. Near the plasma frequency and upper hybrid 
resonance frequency the results are completely useful and a com- 
plicated feature of the experimental result obtained in the 
ionosphere can be explained. This work may be a starting point 
for investigations of the electron acoustic wave effect in wider 
frequency ranges including the cyclotron harmonic effects and 
ion sheath effect. 
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DETERMINATION OF ELF/VLF TRANSMITTING 
ANTENNA PERFORMANCE IN THE IONOSPHERE 

J. Plumer Leiphart 

Satellite Communication Branch 
Communications Sciences Division 
Naval Research Laboratory, Washington, D.C. 20390 

An experiment is planned that will determine the 
effects of the magnetoionic medium on propagation 
between an ELF/VLF source in the ionosphere and a 
surface-based receiver. Proper implementation of 
this propagation experiment requires that instru- 
mentation to measure the driving point impedance 
and efficiency of transmitting electric dipoles and 
loop antennas in the anisotropic ionosphere. 

An ELF/VLF propagation experiment is being planned by the 
Navy to determine how energy emitted from a spacecraft in the 
ionosphere propagates through the ionosphere, couples across the 
air-ionosphere boundary, energizes the Earth-air-ionospheric 
cavity, and propagates to the surface of the Earth. A high power 
transmitter in an omiting spacecraft will be the source of power. 
Receiving and analysis equipment at many geographic locations 
will determine signal strength, polarization and direction of arrival. 
Total time delay and frequency shift will be determined by com- 
parison with telemetered replicas of the transmitted ELF/VLF signals . 

However, losses and other characteristics of the propagation 
path cannot be determined from the received signals unless the 
characteristics of the emitted signals are known. Some information 
relative to the radiated signals can be deduced from measurements 
made on the transmitting antenna. These will be outlined in this 
paper. In addition, the antenna measurements will provide a base 
for future experimentation and guide the design of efficient couplers 
between the high power transmitting source and the terminals of an 
antenna immersed in the ionosphere. 

In order to obtain a wide range of propagation path parameters 
and many "in-situ" conditions for measuring antenna performance, 
a polar orbit with an apogee of 2500 kilometers and a perigee of 350 
kilometers is planned . Including diurnal variations , this will pro- 
vide a range of electron densities from 10 6 electrons/meter 3 to 



■ 
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10 electrons/meter . The predominate ion will change from H to 
+ between the apogee and perigee. The other magnetoicnic para- 
meter affecting the antenna will be the strength of the Earth's mag- 
netic field which will produce electron gyro frequencies of 350 kHz 
to 1500 kHz. 




FIGURE 1 . Angles determining 
sheath characteristics 



/ 

It is anticipated that the ion sheath around the dipole elements 
will be the controlling parameter affecting the antenna reactance at 
low drive levels. This ion sheath will be perturbed by the v x B 
potentials induced in the antenna structure . The size and shape of 
the ion sheath will depend greatly upon the magnitude of v and B 
and the angles 6 , e , and £ , shown in Figure 1 , which uniquely 
define the orientation of the antenna relative to B and v. 

A recent analysis by Mr. Dennis Baker of the Naval Research 
Laboratory based on such a perturbated sheath confirms the antenna 
capacity measurements of the LOFTI II A spacecraft experiment 
which carried a simple dipole. The spinning antenna on LOFTI II A 
is depicted in Figure 2 as wheel spokes in the x-z plane. The spin 
axis S is along the y axis and the x-y plane contains the Earth mag- 
netic vector 3. The spacecraft carried a two-axis magnetometer 
measuring components of the magnetic field parallel to the spin 
axis and along the antenna. When the antenna was in the direction 
of the x-axis , the magnitude of the component of the magnetic f ield 
parallel to the antenna is at a maximum. During the course of the 
experiment the angular displacement between the position for 
minimum capacity and the x-vXis varied. Among other results of 
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VLF Ground Based Measurements over Stratified Antarctic Media 
G. E. Webber and 1. C. Peden 
Electrical Engineering Department 
University of Washington, Seattle, Wn. 98105 

ABSTRACT 

The near-field phase and amplitude broadside to a long hori- 
zontal dipole near Byrd Station, Antarctica, are measured at 12.8 
khz. The resulting data, together with an appropriate mathemati- 
cal model suggest reasonable values of depth and bulk complex 
permittivity for the glacial ice, and a conductivity figure for 
the underlying earth. 

INTRODUCTION 

The phase and amplitude characteristics of the near field 
broadside to a long horizontal dipole at very low frequencies 
constitute a potential source of useful information about the 
nature of the underlying media [Biggs, 1968]. This study concerns 
Antarctic glacial terrain, and makes use of data collected re- 
cently in the Byrd Glacier area to draw tentative conclusions 
about the nature of the sub-surface media. The method suggests 
the depth of the ice shelf is in the neighborhood of 2250 meters, 
that the bulk relative complex dielectric constant of the ice may 
be approximated by the figure 6.0 - jlU, and that the conductivity 
of the underlying earth approaches 0.10 mho per meters. All of 
these figures are shown to be reasonable in light of available 
data from other sources. 

EXPERIMENT AND MATHEMATICAL MODEL 

Two electric dipoles 21 and 8 miles long respectively, lie at 
right angles near the top of a thick ice shelf at an Antarctic 
site 12 miles from Byrd Station; they may be operated indepen- 
dently of each other in the 3-30 khz. frequency range. Although 
they are used primarily for polar ionospheric studies, the follow- 
ing auxiliary experiment was conducted during January 1969. With 
the 21 mile dipole excited at 12.8 khz. and the 8 mile dipole un- 
excited, a vehicle equipped with VLF receiving and telemetry 
instrumentation was driven away from the dipole center along a 
broadside path with respect to th*» axis of the 21 mile antenna. 
It sampled the vertical magnetic field over a total distance of 
11 miles, approximately 3/U wavelength, in the general direction 
of Byrd Station. The field amplitude was measured along this path 
at ground level; simultaneously, the signal was telemetered back 
to the station for phase comparison at the transmitter. These 
data are interpreted in terms of theoretical predictions to be 
described. Further' comparison is made with findings recently 
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will be under control to provide for variation from low levels up to 
the power obtainable from a source capable of supplying better than 
1 kilowatt into the antenna coupling unit. This power will be mo- 
dulated to provide a constant envelope signal or various pulse com- 
binations . This provides the experimenter with a wide range of 
parameters for propagation path experimentation, ionospheric 
excitation, and antenna performance determination. 

The spacecraft will contain a broad spectrum of instrumentation 
to measure ..he parameters affecting the experiment; in particular, 
methods of measuring absolute spacecraft attitude, antenna orient- 
ation relative to the Earth's magnetic field, and devices for 
measuring electron density and ion constituency. 

A prime factor to be known about antennas is the far-field 
pattern. For this purpose, the ground receiving system is quite 
limited. Actually, any practical measurement technique will be a 
compromise. However, other experimenters will have orbiting 
spacecraft containing receiving equipment which will provide some 
indication as to the far field In the ioncsphore. Similarly, one or 
more spacecraft launched and separated from the transmitting space- 
craft could provida information. A receiving rocket launched in a 
trajectory cutting through the radiated field would also provide field 
pattern information. In addition, a receiver in the spacecraft acting 
iJxe a top-side sounder would establish confidence that the signal 
vas radiated by receiving energy reflected or refracted back from 
discontinuities or inhomogenuities . Unfortunately, it is impractical 
to thoroughly scan the entire volume to determine the patterns , how 
they change with radial distances, and the relative power densities 
in the several propagating modes. 

Most information relative to antenna performance will be obtained 
at the transmitting spacecraft. 
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First, the current into the antenna measured at B (shown in 
Figure 3) and the voltage measured at A across the antenna termi- 
nals , together with their relative phases , will be telemetered to 
Earth. This data will be correlated with the environmental data and 
experiment parameters. 

Secondly (as shown in Figure 3), the antenna coupling unit and 
matching transformer will be servo-controlled to obtain maximum 
power transfer. The position of the transformer tap and the reactance 
position will be telemetered to monitor the performance of the servo 
system. This information, providing the range and bandwidth cf the 
servo system is adequate to follow the changes , will confirm the 
results obtained by measurement of voltage, current, and phase. 

It is anticipated that there will be non-linear effects observ- 
able by changing the level of the drive power. Harmonics and other 
frequencies, detectable on the ground or in the spacecraft receiver, 
will be generated and related to the antenna performance behavior. 
This non-linearity will be reflected in the waveforms of the antenna 
current and driving voltage. These waveforms will likely be tele- 
metered to Earth but only for special tests since they require large 
bandwidths in the telemetry channels. In order to conserve tele- 
metry bandwidth, these waveforms (as shown in Figure 4) will be 
processed by a harmonic analyzer on the spacecraft and the results 
relayed in digital form to the telemetry stations . 
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FIGURE 4 . 
Method of 
determining 
non-linear 
characteristics 

Potentials will probably be sufficiently large to cause the sheath 
to collapse at one end of the antenna as shown in Figure 5. Simple 
measurement of reactance across the antenna terminals A and B will 
not be sufficient to reveal this asymmetry. Measurement of the 
capacity between each of the antenna terminals and the body of the 
spacecraft, terminal C, will provide some information relating the 
capacity between the sheath and each of the dtpcle elements . 
Additional information relative to the sheath formation will be ob- 
tained by putting direct current biases between the antenna elements 
and the spacecraft and by an electron-gun on the spacecraft. 
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FIGURE 5 . 
Measurement of 
pulsating sheath 
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Figure 5 shows a static representation of the sheath. When the 
antenna is used for transmission, the driving potential across the 
antenna terminals and along the antenna may equal or exceed the 
v x B potentials. In any case, there will be a sheath which oscil- 
lates from one and of the dipole to the opposite end of the dipole at 
the driving frequency. Measurements at the driving frequency will 
not reveal the full dynamic nature of the sheath. In order to study 
this pulsating sheath, these capacities will be measured at fre- 
quencies of the order of ten to a hundred times that of the driving 
frequency. The representation of the circuit by three capacitors in 
the lower part of Figure 5 is grossly simplified, since the capacity 
is distributed and the ionosphere bounding the sheath is not an 
electrical short but has a complex impedance as diagrammatically 
indicated in the upper part of the figure. In order to determine some 
characteristics of the antenna-sheath-plasma circuit, these im- 
pedances will be measured at several frequencies and a sheath 
impedance derived from the data . 

This experiment is still in the formulation stages and should 
be of great interest to those interested in the effect of plasma and 
anisotropic media on antennas. Any comments to make this a more 
definitive and productive experiment would be appreriated. 



-50- 



VLF Ground Based Measurements over Stratified Antarctic Media 
G. E. Webber and I. C. Peden 
Electrical Engineering Department 
University of Washington, Seattle, Wn. 98105 

ABSTRACT 

The near-field phase and amplitude broadside to a long hori- 
zontal dipole near Byrd Station, Antarctica, are measured at 12.8 
khz. The resulting data, together with an appropriate mathemati- 
cal model suggest reasonable values of depth and bulk complex 
permittivity for the glacial ice, and a conductivity figure for 
the underlying earth. 

INTRODUCTION 

The phase and amplitude characteristics of the near field 
broadside to a long horizontal dipole at very low frequencies 
constitute a potential source of useful information about the 
nature of the underlying media [Biggs, 1968]. This study concerns 
Antarctic glacial terrain, and makes use of data collected re- 
cently in the Byrd Glacier area to draw tentative conclusions 
about the nature of the sub-surface media. The method suggests 
the depth of the ice shelf is in the neighborhood of 2250 meters, 
that the bulk relative complex dielectric constant of the ice may 
be approximated by the figure 6.0 - jl4, and that the conductivity 
of the underlying earth approaches 0.10 mho per meters. All of 
these figures are shown to be reasonable in light of available 
data from other sources. 

EXPERIMENT AND MATHEMATICAL MODEL 

Two electric dipoles 21 and 8 miles long respectively, lie at 
right angles near the top of a thick ice shelf at an Antarctic 
site 12 miles from Byrd Station; they may be operated indepen- 
dently of each other in the 3-30 khz. frequency range. Although 
they are used primarily for polar ionospheric studies, the follow- 
ing auxiliary experiment was conducted during January 1969. With 
the 21 mile dipole excited at 12.8 khz. and the 8 mile dipole urs- 
excited, a vehicle equipped with VLF receiving and telemetry 
instrumentation was driven away from the dipole center along a 
broadside path with respect to th* axis of the 21 mile antenna. 
It sampled the vertical magnetic field over a total distance of 
11 miles, approximately 3/U wavelength, in the general direction 
of Byrd Station. The field amplitude was measured along this path 
at ground level; simultaneously, the signal was telemetered back 
to the station for phase comparison at the transmitter. These 
data are interpreted in terms of theoretical predictions to be 
described. Further comparison is made with findings recently 
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obtained by the U.S. Army Cold Regions Research Laboratories re- 
garding the underlying ice structure in connection with a six inch 
core drilled through the glacier at Byrd. 

The relevant mathematical model is illustrated in Figure 1; 
it comprises a homogeneous, isotropic, lossy dielectric ice layer 
covering a semi-infinite region of imperfectly conducting earth. 
The ice layer lies below a semi-infinite region whose properties 
are those of free space; plane boundaries are assumed throughout. 
Assuming a horizontal line source at the air-ice interface, ex- 
pressions may be derived for the horizontal electric field at this 
boundary in terms of distance from the source. These are well 
known, and have been presented before [Wait, 1962; Collin, I960; 
Biggs and Swarm, 1965]. The vertical magnetic field is of course 
simply related to these expressions through Maxwell's equations. 
Digital computer techniques are used to obtain accurate solutions; 
they are cast in a general form to permit application to a variety 
of natural media, of which ice and underlying earth represent a 
special case. 



21 mile dlpole oir refl | 0n (free space) 




earth region 
<VP Mo) 

Figure I. 

It is shown that two basic modes are the primary contributors 
to the total field at the interface, namely a component represent- 
ing propagation in the ice layer and one propagating along the 
interface itself. An additional weak lateral wave contribution 
results from propagation just below the earth- ice interface. A 
brief mathematical discussion summarizes their origins in terms 
of an integral evaluation in the complex plane. Predictions are 
made regarding the phase and amplitude of the electric field at 
the air-ice interface in terms of a range of possible complex 
dielectric constants and layer depths. The best fit to the 
experimental data is found to occur for the combination of param- 
eters indicated above. Figures 2 and 3 illustrate the agreement 
obtained under these conditions. 
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DISCUSSION AND CONCLUSIONS 

The vertical thickness of the glacier at Byrd Station has 
been measured directly by drilling to the bottom [Gow, Ueda and 
Garfield, 19683; the resulting depth was 2164 meters, which com- 
pares well with our value of 2250 meters. Other investigations 
have yielded depth figures varying up to 2570 meters in the 
neighborhood of Old Byrd Station [Thiel, Bentley, Ostenso, and 
Behrendt; 1959], where the ice-rock interface is believed to be 
quite irregular along a path roughly parallel to, and only a few 
miles from the path traveled by our vehicle. 

A conductivity of 0.1 mho per meter is slightly high for the 
underlying rock, but not unreasonable for the mixture of granite, 
ice and water found at the bottom of the drill hole by Gow, Ueda 
and Garfield. Furthermore, our method is relatively insensitive 
to the assumed conductivity of the underlying medium, so long as 
this is greater than 0.001 mho per meter. 

No experimental data are available regarding the complex 
permittivity of the ice as a function of depth at the drill site 
or elsewhere in the Byrd Glacier area, although the density and 
complex dielectric constant of the surface snow have been measured 
at several different temperatures at 10 khz. [Rogers, 1967]; e = 
e Q (2.7 - j2.7) at - 30° C, with a density of approximately OA 

gm/cc is a typical figure; the complex permittivity decreases with 
temperature, and increases with pressure packing until impermeable 
ice results at a depth of approximately 65 meters [Hatherton, 
1965], Assuming that this parameter is essentially independent of 
pressure over the entire vertical depth of the glacier, and pri- 
marily dependent on temperature and the presence of impurities, it 
is possible to infer some bounds on the expected average figure. 
Simple straight-line approximations to the temperature profile 
measured in the drill hole [Gow, et. al.] yield an average tem- 
perature near -20° C. This temperature, in turn, suggests a lower 
bound of 3 - j6.7 for the complex permittivity in terms of the 
figure for pure ice at 12.8 khz. at this temperature [Watt and 
Maxwell, I960]. Impurities in the glacial ice can be expected to 
increase the dielectric and loss factors. A s>+ep toward account- 
ing for their presence is made by assuming them to be comparable 
to or less than those of the Athabasca Glacier in Alberta, Canada, 
for which Watt and Maxwell give some data. Their figure of 
-j21, measured in the surface ice at 0° C, is taken as an upper 
bound since it is known that the quantity of interest decreases 
at lower temperatures. Moreover, the greater isolation of the 
Byrd Glacier from civilization suggests a lower impurity level 
relative to that at the Athabasca Glacier. On these grounds, our 
value of 6 - jit (e r = 6, tan 6 = 2.33) appears justified. 
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It is concluded that the measurement of the phase and ampli- 
tude of the VLF ground wave can provide information about the 
nature of glacial ice and the sub-surface land boundary. It may 
be viewed as a valuable adjunct to other methods currently in use 
to study these characteristics of the Antarctic terrain. 
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RADIATION FROM AN ELEVATED VERTICAL DIPOLE 
ABOVE AN EXTENDED GROUND SCREEN 

W. T. PATTON 

RCA, MOORESTOWN, NEW JERSEY 

A good inductive ground screen used to provide low angle 
space wave radiation from a vertical dipole can support 
the Sommerfeld surface wave since the phase angle of tht> 
effective surface impedance can exceed 45 degrees and may 
approach 90 degrees for a very good screen. The Sommerfeld 
attenuation functior is important to accurate calculation 
of ground screen effects for screens of large extent. 

The ability of an antenna system to establish long range links via the 
ionosphere is determined by the characteristics of the space wave 
radiation pattern at low elevation angles. For a short vertical dipole 
over a homogeneous earth, the pattern is given by 

Jkhsim|> . _ -jkh sin* . . Z, . \ 

e + R V][ e J sin* - 1QJQ (1) 

= costl* where R^ = Jf_ 



sin<J> 4 



Z l<*> (2) 



is the Fresnel reflection coefficient for vertical polarization. ^1^ 



1 /. cosM 1/2 / , o \l/2 ^1. 

~f n't " IF*) ' n = \ £ r "t ) O) 

is the surface impedance of a finitely conducting earth. The space 
wave radiation pattern (1) of a vertical dipole generally shows sharp 
attenuation for elevation angles below 10 to 15 degrees. 

Two methods of improving low angle radiation are 1) increasing vertical 
directivity by elevating the dipole and 2) providing a ground screen 
under the antenna of sufficiently large extent to be effective at low 
elevation angles. Ground screens longer than 100 wavelengths are 
required to control the pattern at angles about 5 degrees elevation. 
Following Wait (1967) with the Sommerfeld attenuation function included 
as in Wait (1963) the space wave radiation pattern of the composite 
system is givctn approximately by 

ka 



G 3 = G, + 



Z 2- Z l 1+R n / ^ - j Jx 2 +H 2 



J nW e 



1 If 2 J H 2 +X 2 6 vJ(X costy ) W (X.ZJdX 

o i i. 

The Sommerfeld attenuation function is (5) 
W(X,Z 2 ) = 1 -jJTp" e" W erfcU^.w =^1 + | 2 P , P = -j |^ 2 
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and is the composite surface impedance formed by the parallel 
combination of the surface impedance of the earth (equation 3) and the 
impedance of the ground screen. If the ground screen is composed of 
parallel wires or a wire grid, its impedance is given approximately by 



where d is the grid spacing and c is the wire radius. The grid 
impedance is inductive and, to be effective in improving low angle 
radiation, must be much smaller than surface impedance of the ground. 
Thus the composite surface impedance determining the behavior of the 
Sommerfeld attenuation function will be inductive with a phase angle 
approaching 90 degrees. The magnitude and phase of this function shown 
in Figures 1 and 2 show the effect of the Sommerfeld surface wave as 
noticed by Wait (1957) for highly inductive ground systems (phase angle 
over 45 degrees). It is manifest by the attenuation function exceeding 
unit magnitude and undergoing rapid phase variation. 

Satisfactory results can be obtained from equation (4) by setting W 
equal to one only when the magnitude of the numeric distance P in 
equation (5) is much less than one. This may often be the case for 
antenna systems where good radiation below elevation angles of about 
10 to 15 degrees is^of primary concern. However, for good performance 
at lower elevation angles, an extended ground screen must be used and, 
in the interest of economy, the largest wire spacing consistent with 
good performance is desirable. In this case, the Sommerfeld attenuation^ 
function must be included in equation (4). 

The effect of grid spacing on the elevation pattern can be determined 
initially by considering the composite ground system to be of infinite 
extent. The patterns shown in Figure 3 are for this case. The 
interesting result is that for large grid spacings, the radiation 
pattern of the elevated dipole is poorer than that over unaided ground. 
The reason for this can be seen when equation (1) is approximated for 
very small elevation angles 

cosifr sinj / \ 

where the subscript 2 denotes the impedance of the earth as augmented 
by the wire grid. Since is highly inductive, is reduced as the 
dipole height h is increased and obtains a minimum at a height depending 
upon the magnitude and phase of Z-. The pattern of a finite extended 
ground screen will tend to be better than those shown in Figure 3 at the 
lower elevation angles due to diffraction from the end of the screen. 

The Sommerfeld attenuation constant should be included in calculations 
leading to the design and evaluation of economical ground systems for 
low angle radiation. Elevating the dipole above the ground screen can 
in some cases result in poorer performance than might be had with a 
monopole on the ground. 




(6) 
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NEAR FIEI JD SOLUTIONS FOR A VERTICAL HERTZIAN DIPOLE 
OVER A FINITELY CONDUCTING EARTH 

David C. Chang, University of Colorado 
James R. Wait, ESSA Research Laboratories, Boulder, Colorado SO 3 02 



Various approximations to the exact formal integral 
representation are considered. A new approximate form 
is obtained which appears to be valid even when the dipole 
is near a poorly conducting earth. Numerical compari- 
sons with the exact integral formula verify some of the 
conjectures. 

1. INTRODUCTION 

The complete electromagnetic field of a vertical Hertzian 
electric dipole of current moment I Q d s, located at a distance z 
above a homogeneous half -space of conductivity a and dielectric 
constant e can be determined from a Hertz vector II which has 
oaiy a z component [Banos, 1966]. The formal exact solution for 
the time -harmonic problem is expressible in the form 

E-VxVxnT ; H =(o-i«e)Vxn a ; (1) 
z z z z 

n z (p.z) = (- lC 4^° d8 ) K(z, Z Q ;p) ; K = K (p) + K (s) , (Z) 

where the primary influence is 

ikjtz-z^ + p 2 ]* 
K (p, (z,z;p) = ^ _ , {3) 

k [M)^P ! l f 



o o 



. * P 2 1* 

o 

and the secondary influence is 



. . n ik (z+z )[l-o 2 ] . 

K (8) (z, z ; p) = i R(«) e ° ° J <k p c) ; (4) 

O v O C f ,2 1? 

O I* ~ C J 

R(«) = [a 2 (1-c 2 )* - (n 2 - o 2 > fyfh 2 (1-c 2 )* + (n 2 - a 2 >* ] . (5) 
Here, n = [(c - i to e)/(-i« € )]^ is the (comclex) refractive index 



of the earth; £ = 120 tr, the characteristic impedance of the free- 
space; k , the* free -space wave number and a z is the unit vector 
along the z direction. In (l)-{5), the time -dependence e~ i(,,t is 
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assun ed. Since the secondary influence K 8 , in (4), is expressed 
onlv in terms of an integration in the Fourier -transform plane a , 
various asymptotic solutions have been developed for different ranges 
of observation [fianos, 1966]. Specifically, the knowledge of the field 
near the vertical axis (i. e. , p 0) is important in the determination 
of the self -impedance of the source dipole using the e. m. f. method 
(Wait, 1962] as well as in the finding of a better approximate solution 
for the current of a finite linear antenna using the integral equation 
approach [Chang and Wait, 1969]. In the past, appr .ocimate formulas 
which give the expression in closed forms have been developed for 
cases when the source dipole is located far away from the earth 
(k Q (z + z ) t 10) [Jasik, 1961] , and when the earth is highly conductive 
[Wait, 1962; Lavrov and Knyazev, 1965], However, h a more general 
situation, the source and the observation point are near the ground 
which, itself, may not be highly conductive. In this case, the avail- 
able formulas fail to give an adequate description of the physical 
problem. Motivation for a more general analysis has arisen in 
connection with the design of HF antennas for special applications and 
more recently in the remote sensing of the moon's surface. 

2. APPROXIMATE SOLUTION NEAR THE VERTICAL AXIS 

In the following derivation, we make the assumptions: 
p«(z + z Q ) and k Q p« 1 , (6); |n| 2 s 10 « (7) 
Branch cuts at a = * 1 and a = ±n, t in the expressions (4) and 
(5) are so chosen that jarg {a 2 - n 2 ) 5 js vfz and jarg (« 2 - 1) * J < ff/2. 
We now let z + z Q = d and split the integration path in (4) into two 
parts: 



1 * i 



«•(»), t •Z' C a f ~\ o « » ik^d(l-c 2 ) _ . cdc 
K x '(z,z ;p) = i( \+ \ )R(a)e o J (k p e) j 

= i \ [n 2 « - (n 2 - 1 + cV ]/[n 2 a - (n 2 - 1 + e 2 ) 2 ] p [1 - c 2 ] 2 



O a, 

i" 



x e k o da d« + ftin 2 e-(n 2 -l-tt 2 ) 2 ]/[in 2 c + <n 2 -l-c 2 ) 2 ] 
o A 

x e lk o dc J o (k Q p[14 a 2 ) 2 ) do . (8) 

j. 1. 
In deriving (8), a change of variables (1 - c 2 ) 2 - c and (o 2 - l) a -» a 
has been applied to the first and second integrals, respectively. A 
In the first integration, o s 1, we can approximate (n 2 - 1 + c 2 ) 2 
by n, when it is compared with n 2 e , based on (7). This 
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approximation apparently is better than c^|nj J . On the other 
hand, the second integrand contributes significantly only when 
o £ l/(k d) because of the exponential decay term in the integrand. 
Again, based upon (7), we can approximate (n 2 - 1 - e 2 )^ by n . 
Also, because of (6), the Bessel function J (k p[l?a 2 ]£) can be 
approximated by 1. To determine the error involved, we can esti- 
mate the remainder or "error" term: 



a(n , d) = fr in 2 a -<n 2 -l.a 2 )^ _ ino^l n ^o " p[1 + e * 3 *) d « 

~ L in 2 a + (n 2 -l-« 2 }* lna + lJ 
o 

_ f r in 2 6 -(n 2 -l.a 2 )* _ ino-1 l d< 

J L. 2„., 2 , „2\2 ine + lj 

o in a Tin -i-U| 

It ie not difficult to show that this is of the order of l/|nj 4 or 
l/(Jn| a k 2 d 2 ) . Consequently, we can approximate (8) as 

K (8, <z. z q ; p) ss K c (z,z q ) 



k cd 



OB 

pi iok d .-, 

=s i \ [(no-lVW + i)] e ° do + \ [(in c -l)/(in a*l)]e 



-ck d 
o 



This can be recombined to give 



, ik d(l-e 2 ) ? . 



K <z, z ) = i \ {I - zAn(l -a 2 )* + 1]} e ° --^r 

c ° - <l-« 2 )* 

o 

= K (z, z ) + K (z f z ) , (9) 
O Cz o 

P ik d(l-c') z , o* o 

K (z,z) =i\e ° -Z-^h = \ l9 + 9 s • (W) 



and 



r 4 ik d(l-e 2 ) 

K (z.z) = -2i \{l/[n(l-« 2 )*+l]}e ° 

o d-o 2 ) 2 

/ i °* 

,. -ik d/n ik de' . . 

Zi o [ o da* 

= e \ e — r- 

n J a* 

1 +l/n 



i 



) w jhui 1 1" nwm in | t--"i-' , J4a.tuuLu.i,wji»iiiJui"jimwun 
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-ik d/n 

= ~ e ° E,(-ik r) ; r (z + z ) {1 + l/n) . (11) 
n * o o 

where Ej is the exponential integral of the first order [Abramowitz 
and Stegun. 19641 , Substitution of (10) and (11) into (9) yields 

^■♦■o* -ikd/n 
o o 

where the complex distance r is defined in (11). It is interesting to 
note that when jnj » 1 and k^d/Jnj « 1, (12) becomes 

ik (z+z ) 



which agrees with the approximate solution obtained by Wait [1962] 
and Lavrov and Knyazev [19651. 

In the limiting case when jnj — * , a perfectly conducting 
ground, only the first term in (12) remains which corresponds ex- 
actly to the contribution from an image source below the ground. 
On the other hand, if k (z + z ) » 1 , using the asymptotic ex- 
pansion [Abramowitz ami Stegun, 1964] of Ej (x) , we have 

ik (z+z ) ik (z + z ) 

CO _ o o 

K I \ - e 2 e 

c***V ~ k (z + z ) " n k (z+z )(l + l/n) 
o o o o 

ik (z + z ) 



o o 



o o 

which agrees with the approximate solution used by Jasik [ 1961] and 
Berry and Chrisman [1966]. It corresponds to the reflected portion 
of a normally incident wave and is the first term of an asymptotic series 
representation intended for use at higher frequencies [Vfe.it, 1962] . 

3. DISCUSSION 
In order to discus e the various approximate solutions, the ex- 
act solution K^ 8 ), in (4), which is expressed only in terms of an 
integration in a complex transform plane a, is evaluated numerically 
and is compared with our approximate solution K in (12); the one 
used by Lavrov and Knyazev [1965] , K L in (13); and the one used by 
Jasik [1961] and Berry and Chrisman [1966] , K in (14). The earth 
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ground is chosen to have a conductivity of 10 mho/m and a relative 
dielectric constant of 10. The operating frequency is set at 100 MHz, 
giving a refractive index of the earth n = 3.175 + iO. 28 which is not 
particularly large. Both the amplitude and phase of these solutions 
are plotted against the normalized distance k (z + z ), i. e. , the total 
distance of the source dipole and the observation po°nt, in Figures 1 
and 2. When the source and the observation point are close to the 
ground (i.e., k Q (z + z ) s 1.0), both K £ and agree reasonably 
well with the exact solution while Kj doe 8 not. Perhaps 

better than what we assume in the derivation, the error involved is 
less than 10$ so long as k (z + z)^ 0.25, or the total distance is 
longer than 0. 04 of a free -space wavelength! Farther away from 
the ground (i. e. , k (z + z ) & 1), the error in K c is almost negligible. 
In the same range, however, the error in the approximate solution 
K l tends to be worse, while the one with Kj is better. This is 
due to the fact that the error in K j_, is on the order of l/|n z | . 

For all practical purposes, K c can be used adequately for a 
medium with arbitrary refractive index. If, on the other hand, |n| 
is large compared with one, K L is preferred, because of its 
simplicity for application to impedance calculations of low antennas 
over a well conducting earth [Wait, 1969]. 
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PUISE PROPAGATION INTO THE GROUND 



by 

Walter L. Curtis 
The Boeing Company 
Seattle, Wash. 98124 

Abstract 

The transmission of a delta function plane wave into a plane earth 
is determined. The coaplex dielectric properties of the ground 
are modeled using an empirical relationship that exhibits a typi- 
cal variation of conductivity with frequency. Numerical results 
are presented for both high and low loss grounds as a function of 
depth. 

Introduction 

The subject of pulse propagation into the ground has continued to be of 
general interest because of the need to analyse the effects on electron- 
ic systems of high amplitude electromagnetic transients. Such transients 
might be produced by lightning or a nuclear burst. The purpose of this 
paper is to gain sor.e^insight into the shielding properties of the earth 
against pulses. Welt has considered this problem in considerable de- 
tail and gives an excellent suranary of previous work by others. More 
recently Kins an" Harrison^ investigated this problem using a Gaussian- 
shaped pulse. The problem considered here is that of a delta function 
plane wave normally incident on a plane earth. In this way, the re- 
sulting pulse at different depths will exhibit features due to the lossy 
earth which are not masked by frequency limitations of the source pulse. 
In addition, the results for any other incident pulF^ shape is just the 
convolution of the incident pulse and the result gi./en here. Previous 
investigations by others have assumed the earth to have a conductivity 
that is constant for all frequencies. Trie present analysis uses a model 
fcr conductivity that more closely matches experimental measurements. 

Analysis 

The problem is solved by finding the steady state solution and then 
taking the inverse transform numerically on a computer to find the 
pulse shape at different depths in the ground. As shown in Figure 1, 
the ' trait emitted wave is in the earth. It is assumed that the 

earth has a permeability of free space and a complex relative dielectric 
constant, E" 5 , given by 

€ ? ~ € i ~ J cj €. B 

where £ 9 , is the usunl relative dielectric constant which as a value 
around 8 to 25 depending on the moisture content, a~ is the conductivi- 
ty of the ground, u> is 2v times frequency, and is the permittivity of 
free space. Noraally the conductivity is assumed constant with the 
frequency, but measurements on typical soil show considerable variation 
with frequency— especially at higher frequencies where it tends to be 
directly proportional to frequency. This is well known for low loss di- 
electrics where the loss tangent defined as oyio^e^ is usually assumed 



independent of frequency. To account for this, the followir.g model for 
T 9 was used in this study: 

where <r$ is the usual low frequency conductivity and typically has 
values from .002 to .02 tahos per mf»ter and L ? is the high frequency 
loss tangent. There is little information available on values of 
for different soils but measurements on a soil sample from the Seattle 
area show a value of 0.5. The above model is good for all practical 
purposes but is not too convenient when extended to the entire complex 
frequency domain. An alternative model that overcomes this and still 
has about the same variation with real frequencies is 



The steady state solution for the electric fiela in the earth, E t (w), 
in terms of the incident field, E t (u) , is 

E t = T e rd E; 

where T is the transmission coefficient of the earth-air interface, Y 
is the propagation constant of the earth and d is the depth, T and 
are well known and can be- written in terms of the complex dielectric 
constant as 

T = " * 

and 



Assuming the incident pul3e in time is an unit amplitude delta function, 
then E. (w) = 1. The field in the earth as a function of time and depth 
is then the inverse transform of E, , i.e., 



e 



VV> = tw) £ i e du > = ^"J i + fr, 

Results -«> 

The above transform has been evaluated numerically on a computer. The 
results for two different conductivities and several different depths 
are shown in Figures 2 and 3 * Figure 2 shows the results for <^ - .002. 
£ ? •=■ 10 and L 9 = • 5 " which is representative of a moist low loss 
soil found in the Seattle area. At all depths, the" pulse is character- 
ized by a positive high amplitude narrow pulse followed by a low ampli- 
tude slowly decaying negative tail. The peak of the pulse shows a 
delay corresponding to that expected for a wave propagating to the same 
depth in loss less dielectric. It should also be noted that the posi- 
tive pulse appears almost symmetrical about its peak. Figure 3 shows 
the results for ooiy a change in the low frequency conductivity, i.e., 
<7£ = .02, <»j = 10 and L 5 = .5* This is representative of a soil 
with about the same moisture content as before, but with a greater 
amount of ionizing material present. It is not known how L 3 might 
change sw it was left fixed for this example , The same pulse character- 
istics are observed but with much lower amplitude. At depths less than 



2 meters, the positive pulse is still quite symmetrical and has a time 
delay the 3ame as before. However, at greater depths the pulse becomes 
skewed and its peak occurs at a greater time delay. For example, at a 
depth of 10 meters the peak occurs at abnut 250 nanoseconds instead of 
100 as before. In addition, the peak amplitude decays much more rapidly 
with depth than it did with the lower conductivity. 

Conclusions 

The impulse response for propagation into the ground has been found and 
typical examples shown. These results show the limiting values of pulse 
width and rise time for any source. For example, at a depth of 1 meter 
the maximum rise time and minimum pulse width will be on the order of a 
ftsv nanoseconds while at 10 meters depth it will be on the order of 50 
to 100 nanoseconds. The results for any other incident pulse shape can 
be found by convolution with the results presented here. 
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PERFORMANCE OF VALLEY SPANNING VLF ANTENNAS 
BASED ON ELECTROMAGNETIC MODELLING 



A. N. Smith and E. J. Jackson 
Westinghouse Georesearch Laboratories 
Boulder, Colorado 



Passive electrical characteristics of numerous transmitting 
antenna configurations, similar to valley- spanning proto- 
types of interest for moderate power application in the 
OMEGA Navigation System, were measured for models on 
alterable valleys. Resulting parametric curves and 
modelled near-zone field distributions plus soil con- 
ductivities of candidate sites enable accurate prediction of 
power and bandwidth capability for cost studies. 

The communicator/navigator generally hands his problem to the 
communications systems engineei in the form of a desired reliability 
(availability in some noise environment) in a coverage area for some 
signal format yielding a specified information rate. These require- 
ments filter down to the antenna system designer as a minimum 
radiated power (P y ) at some frequency (f) or range thereof in some 
specified minimum bandwidth (B ). Present-day technology supplies 
an upper limit on operating ✓oltage, V., and previous generalized 
studies give a range of desirable efficiency r\ (and hence input power 
P ) functionally related to radiated power level so that costs will 
probably be near a minimum. For the OMEGA System, the respec- 
tive values have been established as P >: 10 kW, 10. 2 < f < 13. 6 
kHz, B > 10 Hz @ 10. 2 kHz; V. < 25% kV, P ^ 150 kW; ti 2: 
6.7%. aS 

At VLF, where transmitting antennas are almost unavoidably electri- 
cally short, these six variables and four others are related by the 
four equations given below (applicable in general to any simple tuned 
circuit of lumped parameters and specifically to the one illustrated): 

P = 6.95X10~ 43 C^ h 2 V 2 f 4 H I 2 R 
r e t a r 

B = 1.108 X 10" 13 C A h 2 ti " 4 f 4 = f/Q 
as e as as 

P = P /ti (MKSQ units) 

v" = [ 1 - F(f/f ) 2 ] V. (f Q = (2. yiTc^)- 1 ) 



The remaining four quantities C^ (electrostatic capacity), f (self- 
resonant frequency), Vj. (mean topload voltage) and h g (effective 
height) are thus determined if the original six are fixed; this is 
usually not the case, but certain ones are allowed to take on values 
below well -defined upper limits, so that cost- and performance -trade 
comparisons can be carried out for several candidate sites as well as 
for several configurations in final single selection. A significant 
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possibilities for tuning helix resistance against ground system 
resistance R , which are the two major loss components, always in 
regard to the^condition of antenna system efficiency. In performance 
trading, one finds besides that specification of both power and band- 
width at more than one frequency in the operating band over- 
determines the system. Thus one specifies minimum acceptable for 
each at that frequency it is hardest to attain; e.g. , for OMEGA 
P at 10. 2 kHz, and B (under an earlier concept) at 13. 6, but not 
both at both frequencies. V must be corona free. 




h / 




t? ■ R/R 

as r os 



a: antenna 
as: antenna system 



R « 1607T 2 (h A) 2 ; 0.02^R ^0.10 
r e r 



OT R h -V Q h' V"< X a- 1/wC 0> 



R +R d £0.07a f° ; R h a f" 1 ; R^f 2 



Matching 



W WW 

1 R R, 
r 9 a 



c: conduction Q is nearly 
d: dielectric constant 



g: grouna 
r: radiation 



h: helix + vario. 



The environment is significant to performance determination in two 

respects: geometrical, relating to C^, h {capacity centroid), and 

f (distance relative to wavelength over which charge muslbe carried 

to components of C Q ); and electrical, relating to the nature of the 

impedance plane, losses in which must be controlled by design of 

some ground system so as to simultaneously permit attainment of 

P and B . Since the cost optimization process must consider both 
r sis 

fixed (installed) as well as operating (capitalised over stated period 
of years), some latitude is offered to charge against site modification 
for efficiency (such as elevating spans on towers) versus more trans- 
mitter and prime power. 



The power, bandwidth, voltage, and efficiency figures given above 
imply, through the fundamental equations, that we look for sites 
yielding effective heights for OMEGA applications, from 100 to 200 
meters and capacity from 0.050 to 0.030 ufds; the inverse relation- 
ship of capacity and inductance usually results in f Q between 25 and 
40 kHz, so that V is within 10% of V. at the low end of the band. 
Reactances to tune are then such as to require helix Q*s between 1500 
and 3000 in order to hold R^ to values such that the residual ground 

system loss R can be obtained with reasonable installation. 
y g 

The valley -spanning approach to obtain these figures is advantageous 
only because of the relatively modest power and efficiency necessary in 
the OMEGA system, since the very presence of the ridges uced as 



supports for the spans precludes economic installation of an 
extremely elaborate ground system and also results in high field 
(both H- aDd E- fields) concentrations in locations difficult to reach 
with a wire grid. A high power, high efficiency installation requires 
the span supports to be such as not to get in the way. Displaying the 
loss budget computation for the ground system iu the form 

R = (R TT + R_). . , + (R T + R_) , + R . + R ... 
g H E inside H E outside wires termination 

in which 

R *f R_ -- Re (l/I 2 ) \ [t, (HJ 2 + — w 2 « 2 E 2 ] dA inside: a. to a 
HE a^JOcuz io 

outside: a to 1/6 „ 
o 

shows the critical role that field distributions play in making perform- 
ance predictions have to do with loss budget, The most convenient 
way to arrive at these is through scale modelling. 

Capacity can be calculated rather easily to a fair degree of accuracy, 
because its logarithmic dependence on the ratio h/a makes this rather 
insensitive to inaccuracies in knowledge of height; not so with effective 
height, since it is a summation of production of local height with 
capacity. Thus again modelling by analogue methods is of great aid 
in accurate performance predictions. 

Some 60 configurations for valley-spanning antennas for the OMEGA 
application have been modelled at WGL, at 480/1 scale in air. 26 
more were done for one of the high-power communication stations. 
Another 1100 variants were studied in an electrolytic tank. The 
results generally show that the required C n and h p for OMEGA power 
and efficiency can be obtained with from 6 to 4 spans over valleys 
7, 000 or more feet wide, with aspect ratio 4/1, without use of 
support towers, and with ground, systems of the order of 2,000 foot 
radius if the conductivity is 10" J mhos per meter or higher. 
Extremely high soil conductivity is not advantageous if there is a. 
large seasonal decrease in effective conductivity {e.g. snow) for 
which provision must be made in ground system design to meet 
radiated power requirement because danger exists that during high- 
conductivity portion of the year bandwidth cannot be met without 
artifically supplying extra loss. From field strength pattern meas- 
urements, no evidence exists that valley-spaa, antennas of the 
proportions studied exhibit any characteristics of slots, as was one 
time discussed. 



Work performed under sponsorship of OMEGA Project Office, PM-9, 
through Chesapeake Division, Bureau of Yards & Docks, U. S. Navy, 
under Contract N62477 - 68-C-0940. 
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SURFACE WAVE CONTRIBUTIONS TO THE BACKSCATTERING 
FROM LARGE, HIGH- DENSITY DIELECTRIC SPHERES. 



BY 

M. A. PLONUS AND H. INADA 
DEPARTMENT OF ELECTRICAL ENGINEERING 
NORTHWESTERN UNIVERSITY 
EVANSTON , ILLINOIS 



An application of the Watson transformation to the Mie series 
divides this solution in geometric optics and diffracted fields. 
It will be shown that the geometric optics fields are smaller than 
heretofore assumed. The diffracted field comes from two types of 
surface waves. One of these is the dominant contributor to the 
backscatter. Numerical results for values of ka from 5 to 20 will 
be presented. 



The following figure shows the backscattering from a dielectric sphere 
as calculated from the Mie series. The top continuous curve shows the re- 
sults which are obtained when the geometrical optics method is applied to 
the sphere [Atlas and Glover, 1962]. There is little agreement with the 
exact results. After an application of the Watson transformation to the 
Mie series an integral can be split off which can be identified with the 
geometric optics contribution since it comes from the lit region of the 
sphere [ir.ada and Plonus, 1969]. The lower curve in Fig. 1 shows the re- 
sults which were obtained from this contribution. The geometric optics 
contribution for this particular range of ka clearly cannot account for the 
major backscatter. 




Fig. 1. Backscattering 
Cross section curves 
from the Mie Series and 
from Geometric Optics. 
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Tbe Watson transformation yields other terms which are residue contribu- 
tions and can be identified with surface waves. These waves can travel about 
the sphere many times. Fig. 2 shows the residue contribution of one pole 
(mode 1). N=l is a wave that has travelled around the back of the sphere 
once, N-2 is a wave that has travelled twice, etc. 




Fig. 2. Surface Wave 
Contribution to Back- 
scattering (dominant 
pole only) . 



Once the existence of the surface waves has been established, it is not 
difficult to see that interference between waves which travelled around once, 
twice, etc. can exist. 

ihe next figure shows the interference curves for all N for three separate 
poles (1-1,2,3). Let us examine the resultant curve due to the first pole. Up 
to ka=7 we have a simple periodic curve which when we re-examine Fig. 2 could be 
the interference between N=l and N=2. For ka> 7, the 1=1 curve becomes doubly 
periodic which means N=3 is coming into play. Similar observations can be made 
for the other pole contributions. Notice that for higher l's the corresponding 
residue contribution has a smaller peak which occurs at larger ka. 




Fig. 3 Surface Wave 
Contribution to Back- 
scattering showing inter- 
ference between surface 
waves of N=l,2.3- 
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The next figure shows that the surface waves are the major contributors 
to the backscattering. The small periodicity of ka » .7 can be explained as 
being due to interference between surface waves which have gone around the 
sphere a number of times, whereas the large periodicity of ka « 10 seems to 
be due to the contribution of the different poles (1=1,2, •••) which peak at 
different values of ka. 




Fig. 4. Backscattering 
Curves showing that the 
surface waves are the 
dominant contributors 
to the total backscatter. 



H. Inada and M. A. Plonus (1969), "Numerical Results for Backscattering from 
Large, High-Density Dielectric Spheres", Proc. IEEE, 51, No. 6, 1192-1193. 

D. Atlas and K. M. Glover (1962), "Backscatter by Dielectric Spheres with and 
without Metal Caps", Proc. Inter-Disciplinary Conf. Electromagnetic Scattering 
(Pergamon Press, New York), 213-236. 
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MICROWAVE TRANSMISSION MEASUREMENTS 
THROUGH RADIANTLY HEATED ANTENNA WINDOWS 

A. C. Lind, J. F. Fox 
McDonnell Douglas Astronautics Company - Eastern Division 

and 

D. Q. Durant 
McDonnell Aircraft Company 
St. Louis, Missouri 



A new technique has been developed which allows continuous microwave 
tranp-nission measurements to be made through reentry vehicle antenna 
winders while they are subjected to high heat fluxes. The antenna 
window specimen is radiantly heated by a cylindrical graphite heating 
element which produces a heat flux of approximately 280 Btu/ft^-sec. 
Unique data is shown on the degrading effect of moisture absorption 
in slip-cast fused silica. 

1. Introduction - During reentry, the antenna window of a missile is sub- 
jected to very high heat fluxes. Communication system performance is degraded, 
not only by the plasma that forms around the vehicle, but also by losses 
induced in the window material as a result of the severe heating. Dielectric 
constant and loss tangent data is generally not available in the temperature 
region of interest, that is, greater than 3,000°F. As a result, others (Mead, 
1967; Cetaruk, 1968) have performed microwave transmission measurements 
through antenna window samples that have been heated with a plasma jet. This 
type of measurement provides a qualitative estimate of antenna window perfor- 
mance during reentry. In the plasma jet test, the transmission measurement is 
possible only after plasma jet shutdown because of the interference effects 

of the plasma during the heating cycle. The interpretation of such results 
can be misleading if transient material effects are important or if the window 
cools significantly in the finite time required to shut down the plasma jet. 
The facility described in this paper employs radiant heat and is designed so 
that the heat.'ng element does not interfere with the microwave measurement. 
Continuous microwave data is obtained during the entire heating and cooling 
cycles . 

2. Description of Facility - A schematic of the facility is shown in Figure 1. 
The antenna window specimen is supported in a water cooled, open-ended X-band 
waveguide that transmits to a receiving antenna located several feet away. The 
waveguide is located at the base of a tubular, electrically heated (U5 kW), 
graphite element that is used to radiantly heat the window specimen. The sys- 
tem is purged with ritrogen to minimize burning of the graphite element which 
is at a temperature of approximately 5>000°F. A loose fit between the specimen 
and waveguide results in poor heat transfer and minimizes cooling of the win low 
by the waveguide. The specimen also protrudes slightly out of the waveguide to 
insure a uniform surface temperature. The tubular heating element concept 
allows continuous microwave transmission during the heating cycle. 



- ( I - 



The waveguide assembly is interchangeable with an identical unit containing 
a water-cooled copper calorimeter to obtain heat flux calibrations. A typical 
heat flux profile is showu in Figure 2. A steady state heat flux near 
280 Btu/ft^-sec is reached after approximately 30 seconds. The test time is 
arbitrary, and for the data shown here was 2 minutes; this test time is suf- 
ficient to produce a significant melt layer on slip-cast fused silica window 
specimens. 

3. Microwave Transmission Data - The received microwave signal is recorded on 
an X-Y recorder which is calibrated with a precision attenuator before each 
test in terms of db loss. Directional couplers are employed in the transmis- 
sion leg to allow visual monitoring of the signal generator output and reflected 
power level during the test. Typical transmission data at a frequency of 10 GHz 
are shown in Figure 3 for slip-cast fused silica. The top trace was recorded 
with no window in the waveguide, and demonstrates the insignificant effect of 
the heating element on the microwave measurement. Two specimens taken fi*om the 
same block of material were tested at the heating rate shown in Figure 2. One 
specimen was dried in a low temperature oven and then stored in a desiccator 
prior to testing. The other specimen contained a normal amount of moisture 
absorption. The dry specimen showed a gradual decrease in signal level reaching 
a steady state attenuation of 0.5 db. The specimen containing the moisture 
showed a sharp decrease in signal level corresponding to an attenuation of 7 db 
and then gradually returned to the attenuation level of the dry specimen. 
These test results demonstrate the importance of continuously monitoring trans- 
mission measurements during high temperature testing. No significant change 
in reflected power level was observed during either test. This is in contrast 
to the results of plasma jet tests where mechanical erosion of the window occurs, 
causing detuning changes that are difficult to take into account during data 
reduction. 

k. Summary - The technique presented here is an improvement over other methods 
since it allows continuous microwave measurements during the entire heating 
cycle. The data presented for slip-cast fused silica demonstrates the impor- 
tance of continuous measurements, and shows the degrading effect of moisture 
absorption during the initial heating of the material, which could affect 
communication system performance in a reentry situation. Also 4 this technique 
does not produce mechanical erosion of the window, so that the measurement of 
window material losses are not masked by detuning effects. 
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MEASURED HEAT FLUX AT ANTENNA WINDOW LOCATION 
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FIGURE 2 



EFFECT OF MOISTURE IN SLIP CAST 
FUSED SILICA WINDOWS 




10 20 30 W SO r IM 130 

Tl«£ . SEC 



FIGURE 3 
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OPTIMUM FREQUENCY FOR A VLF T E LEC OMMU NIC ATJON 
SYSTEM USING BURIED ANTENNAS 



R. Gabillard, J. Fontaine, P. Degauque 
University of Lille, France 



Abstract 



This paper is a study of the best choice of frequency for a VLF 
telecommunication system between buried receiving and transmitting 
antennas. The efficiency of such a system depends on numerous 
parameters. We show that there is an optimal direction of the em- 
bedded antennas which is different according to we use a frequency 
higher or lower than a critical one, . We al3o show that the 
electric field E r increases toward a maximum value at a frequency 
*opt • ^ e §i ve the expressions of f c and fopt . The ratio of the 
electric field Er at the optimal frequency to the electric field E<p 
at ELF brings into evidence the interest of the choice of the optimal 
frequency to establish a VLF subsurface communication. 



PROPAGATION BETWEEN PARttLLEL IMPEDANCE SURFACES 



R. B. Dybdal 

The Aerospace Corporation, El Segundo, California 

The modal structure of a parallel plate waveguide characterized 
by impedance boundary conditions is analyzed. A formal solution 
for both parallel and perpendicular modes may be obtained and 
impedance loci are given for the lossless equal impedance case. 



The propagation between parallel impedance surfaces characterized by 
impedance boundary conditions is analyzed. This problem was first 
considered by Barlow and Cullen [1953] and later by Barlow [1965j and 
Wait [1967]. The geometry is depicted in figure 1. The modal structure 



' ' ' ' * ' ' f - 
^•Zg y = H- 



///////// /*/ 



FIGURE 1. Geometry of Parallel Plate Configuration 

may be formulated from Hertzian potentials [Collin, 1960J and a suitable 
potential function is 

(1) *H,E - (e^ + R H ,E " y) > e" j(l? + 

The subscript H goes with perpendicular polarized modes; the subscript 
E, with parallel polarized modes. Real values of p Q correspond to bound 
surface wave type modes; imaginary values, to plane wave components 
bouncing between the surfaces. 

The solution for the parameters p and Rh e for arbitrary impedance 
values has been obtained and the orthogonality relations reveal modal 
energies are coupled in the lossy case [ Eybdal, 1968) . For equal value 
impedances, Rg and Rh reduce to ±1 simplifying the modal solution. The 
properties of this case are summarized in table 1. Impedance loci may 
be devsloped for the lossless case. Note that losses make p complex. The 
impedance loci are depicted in figures 2 and 3 and p Q is obtained from 
the intersection of the loci and a horizontal line drawn at the surface 
impedance value. The usual impedance requirement for bound surface wave 
modes is demonstrated here; i.e. inductive impedance is required for 
bound perpendicular modes; capacitive, for bound parallel modes. The 
cutoff condition for bounce modes is jp Q { = k and is shown in the loci 
for a waveguide 1.3 wavelengths in height. 
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TABLE 1. Mode Characteristics for the Equal Impedance, Lossless Case 
Mode R EjH Nature of p Q Impedance Value Nature of Impedance 

Perpendicular Modes 

.rMtanh-^S. inductive 
_jJ|altanll^ 
.il2oicothl£aiS- 



A 
B 
c 

D 



+1 
+1 
-1 
-1 



real 

imfu&nary 
real 
imaginary 



«Lcoth» ^ 
we ? 

jkoLcot -LBaiL 

W £ 2 



capacitive- 
inductive 
inductive 



E 
F 
G 
H 



+1 
+1 
-1 
-1 



real 

imaginary 
real 

imaginary 



Parallel Modes 

-jjjLcothlftjH 

, uju .. JPo H 
^{Pot 00 * 2 
-AtanhbgLS. 

-jrr^rtan 



tPnlH 



inductive- 
capacitive 

capacitive 

inductive- 

capacitive 
capacitive 

capacitive- 
inductive 




FIGURE 2. Impedance Loci for Perpendicular Modes 
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FIGURE 3. Impedance Loci for Parallel Modes 



This work was performed, while the author attended The Ohio State 
University and was supported in part by a sub-contract between The 
Andrew Corporation and The Ohio State University Research Foundation, 
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IME3DAITCS OF A DIPOL5 IN A LAMB LABORATORY rAGNSTOPLASMA. 

S.Y,X. Tarn 

Plasma and Space Physics Laboratory 
RCA Limited 
Montreal, Canada 

■ 

Impedance of a 9.2 GHz ^-dipole in a large dc magnetoplasma 
(diameter « 7 ao ) have been measured for various electron densities, 

e<Xsf-2S) * 2, and magnetic fields, 0<ysj~\ < 1.5, Ao being 

V a / \ / 

the frce-spaca wavelength, w, u , and w the si^inal, plasma and 

pe ce 

electron cyclotron frequencies respectively. Results have been 
consistent with the theory. 
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INTEGRAL EQUATION APPROACH TO THE RADIATION FROM A 
VERTICAL MONOPOLE OVER AN INHOMOGE NEOUS GROUND PLANE 
(To appear in January 1970 issue of Radio Science ) 

J. R. Wait and K. P. Spies 
ESSA Research Laboratories, Boulder, Colorado 80302 



The possibility that HF antenna radiation can be enhanced at low 
angles, by a tapered ground system is explored. The model consists 
of a vertical linear antenna erecte 1 over a radially inhomogeneous 
ground system. The surface impedance xb chosen to be a smooth 
monotonic function of the distance from the base of the antenna. In 
a practical scheme, this can be achieved by locating a radial wire 
system in a symmetrical fashion beneath the monopole. The com- 
bination of the equivalent wire grid impedance and the ground im- 
pedance is computable from previously developed theory. 

It is first shown that the tangential magnetic field satisfies an 
integral equation of the Voltsrra type. Allowing the effective surface 
impedance to vary in an exponential fashion, a power series solution 
for the tangential magnetic field is then obtained. ThiB involves 
integral and half-integral powers of the distance from the bass of the 
antenna to a point on the ground surface. As a check for the homo- 
geneous ground (i.e. , no wire grid system), it reduces to the 
appropriate series expansions for the Sommerfeld attenuation function. 
The power series solution, so obtained, is supplemented by a 
numerical evaluation of the integral equation. 

Using the solutions obtained above tor the tangential fields, the 
radiation pattern of the antenna is calculated for several inhomogeneous 
ground models. Some earlier results using simplified assumptions are 
confirmed. + It is shown that a tapered ground screen model is mucl 
better than an abruptly truncated system. Some consideration is also 
given to intentional azimuthal variation of the ground system. 



J. R. Wait (1967), Pattern of a linear antenna erecte? over a tapered 
ground screen, Can. J. Phys. 45 (9), 3091-3101. 
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M1CR0WAVE (4.75GHz) SURFACE IMPEDANCE MEASUREMENTS 
OF A 1* DIELECTRIC WEDGE ON A PERFECT CONDUCTOR 

R.J. King and C.U. H us tig 
University of Wisconsin, Madison, 53706 

Abstract: 

When propagating in the direction of increasing depth, the 
results near the wedge apex compare best with calculations based 
upon simple parallel stratified media theory. When propagating 
in the direction of decreasing wedge depth, the results compare 
reasonably well with the theory given in 1967 by Sehlak and Wait. 
However, the nature of the surface impedance is decidedly dif- 
ferent for the two cases, being inductive near the wedge apex for 
the first case, and capacitive with a negative real part for the 
second case. 
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SESSION VIII, 1:30 p.m., Thursday, July 17, 1969 
Round Table Discussion on: 
DESIGN TECHNIQUES FOR PATTERN CONTROL BY GROUND SCREENS 

SYNOPSIS by R. V. Row, Sylvania, Waltham, Mass. 



The moderator (R. V. Row) introduced the two others sitting 
at the table with him: Dr. J. R. Wait (ITS/ESSA Boulder, 
Colo.), and R. D. Wengenroth (General Electric HMES, Syra- 
cuse, N.Y. ), and proceeded to invite the participants to 
range further afield in their discussions than the adver- 
tised title of the round table would indicate and to con- 
sider terrain, fences and other obstacles within the scope 
of the meeting. 

The discussions were then opened by J. R. Wait who dis- 
coursed on the proper use of transfer impedance and input 
impedance between antennas. He pointed out that most radio 
propagation work in the literature deals only with the 
transfer impedance and when terms like transmission loss or 
system loss are used by an author care must be taken by the 
reader to ascertain how the input power was calculated. 
Frequently input power is calculated on the basis of a di- 
pole in free space. The speaker then proceeded to sketch 
briefly how a perturbation technique, often known as the 
compensation theorem, has been used to compute both input 
impedance and transfer impedance. For input impedance of 
an antenna and ground screen one need only be concerned 
with fields close in to the antenna, whereas the fields over 
the entire ground surface lying between the antennas contri- 
bute to the mutual impedance. He commented on the wide 
utility of the method aad the fact that it has variational 
stationary properties which reduce the relative error in the 
perturbed quantity being calculated, to a lower order than 
thb error in the assumed tangential E or H fields. 

The moderator invited E. W. Seeley (U.S. Naval Weapons Cen- 
ter, Corona, Calif.) to comment on his calculation of the 
radiation efficiency of a long horizontal above ground wire 
antenna. (See pp. 174-178 of Vol. I of EEAP Proceedings.) 
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It was stated by See ley that his expression for radiation 
efficiency is based on the actual input power to the an- 
tenna. The moderator attempted to restate this in terms of 
the E field maintained at a point in space by the actual 
antenna and a reference dipole each with the same input 
power. Seeley stated that it is more useful at VLF to con- 
sider power density radiated in a certain limited range of 
angles (10-15 degrees) above the horizon, rather than total 
power radiated, since this is the power that maintains the 
field at a distance via ionospheric reflection. 

The moderator invited W. L. Curtis (Boeing Co., Seattle, 
Wash. ) to comment on the experimental and interpretive 
theoretical work done by he and Coe [ 1964] on antennas on 
large ground screens over lossy earth. He pointed out that 
both the experimental fork and a subsequent theoretical 
study [ Curtis 1964} using the compensation theorem showed 
that for engineering purposes the radiation pattern of any 
low profile antenna on a finite size ground plane on a flat 
earth can be easily obtained from a free-space pattern of 
the antenna and ground plane. This means that modeling or 
other free-space antenna design techniques can be used, even 
if the antenna is to be placed over flat ground. To first 
order one earth loss correction curve can be used for all 
frequencies above the HF band and all ground planes greater 
than A/2 in radius. The small change in this 'universal' 
correction due to ground screen size can be taken into ac- 
count using the curveu in Curtis [ 1964] . In the HF band 
and at lower frequencies where the ground refractive index 
is frequency dependent different correction curves would 
need to be prepared. 

Professor L. B. Felsen (Brooklyn Polytechnic Institute) then 
suggested that there is another way of looking at antennas 
and ground screens, namely ray optics which provides useful 
insight not only jfc c the input impedance but also for the 
radiation field. He illustrated his remarks by the example 
of a radiator over the center of a circular screen. In 
this case all the edge scattered rays travel in radial dir- 
ections, and the scattering properties of an edge are well 
known and available in the literature. This way of looking 
at the problem also provides insight into the effect of 
non-circular screens, for example an elliptical screen, 
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where only four rays are scattered back to the center. 

T. Kaliszewski (General Electric HMED, Syracuse, N.Y. ) 
pointed out that diffraction and ray optics concepts were 
combined by Senior [ 1956] in studying ground wave propaga- 
tion across edges. Krause [ 1967) has developed similar 
ideas for looking at sloping perfectly conducting screens 
and non-planar screens. Apparently the method is not 
readily adapted to imperfectly conducting screens, which is 
the situation of greatest practical interest. 

Professor Felsen suggested the the surface impedance con- 
cept could be used with imperfectly conducting screens. 

V. Arens (Sylvania Electronic Systems, Mountain View, Calif.) 
commented that he had adopted this point of view and used the 
plane wave reflection coefficient appropriate to the specular 
point in question to calculate the radiation pattern from a 
finite vertical dipole over an imperfect ground screen. 
Edge diffraction was not considered and is believed to be a 
second order effect [ reported at the ARPA sponsored OHD sym- 
posium, Boulder, Coir 1966] . Comparison was excellent 
between predictions m ,e this way and measurements made on 
a large circular array in the HF band. 

J. R. Wait stated that indeed the ray approach gives in- 
sight into the pattern effects of ground screens. He pro- 
ceeded to explain an approach he used to study the pattern 
of an array of slot antennas in a semi-infinite ground plane 
above a lossy half space. Employing the reciprocity principle 
he considered the incident Wave (vertically polarized) from 
above, removed the lower dielectric half space and replaced 
its effect by an image plane wave incident from below. This 
procedure satisfies the boundary conditions except in a 
region close to the edge of the half-screen along the die- 
lectric interface and indeed the results agree with those 
found by Clemmow [ 1953] who used a somewhat similar plane 
wave spectrum approach. The results also agree with those 
found from the compensation theorem. [These results were 
presented at the Congres International Circuits et Antennes 
Hyperfrequences, Paris, October 1957. t Experimental results 
dealing with this same problem were reported at this same 

t ~ * , 

J. R. Wait, L' Onde Electrique, 38e Annee, No. 376 bis, Supplement 
Special, vol. 1, pp. 21-29, 1957. 
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congress by Robieux and Simon.] The speaker then went on 
to formulate an exact two dimensional Integral equation 
for inhomogeneous ground problems and cast in the form of 
attenuation functions. He also suggested that the ray 
optics approach is not easy to justify in many practical 
cases where only modest or even small size ground screens 
are being considered. 

W. L. Curtis commented that the ray theory shows that for 
a truncated ground plane the radiation pattern has small 
oscillatory lobes just where measurement and the compensa- 
tion theory say they are, and that this kind of 'prediction' 
works even on planes as small as a half wavelength in 
radius . 

The moderator then invited E. A. Thowless (U.S. Naval 
Electronics Laboratory Center, San Diego, Calif. ) to re- 
view the work they have done on ground screens. Thowless 
briefly reviewed the contents of four NEL reports as 
follows: 

NEL Report 1346 - "Ground System Effect on High-frequency 
Antenna Propagation," by W. E. Gustavson, W. M. Chas? and 
N. H. Balli, 4 January 1966. Based on model measurements 
it shows that there is significant improvement (5 to 7 dB) 
in radiation pattern at angles below 20 degrees due to 
ground screens of radius 2 1/2 to 7 1/2 wavelengths and 
moderately conducting ground at frequencies of 10 to 30 MHz. 
It was deduced that mesh sizes of 1 or 2 feet would be ade- 
quate for screens used in the HF band. 

NEL Report 1359 - "HF Extended Ground Systems: Results of 
a Numerical Analysis," by G. D. Bernard, W. E. Gustavson 
and W. M. Chase, 24 February 1966. Based on theoretical 
work of Wait and Walters [ 1963] on circular sector ground 
screens. Computations are presented concerning circular 
screens for frequencies of 4, 8, 16, 32 MHz, ground screen 
radii of 2, 4, 8, 16, 32, 64, 128 wavelengths, grid spacing 
of 6, 12, 24, 48 inches and ground parameters of € = 4, 
a = .003 mho/m. € = 10, a = 0.01, € = 20, a = .03 and € = 30, 
a = .1 and elevation angles of 2, 5, 10, 15, 20, 25 degrees. 
In general an improvement of 3 to 9 dB is attainable over 
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poor soil with a ground screen of 250 foot radius. The 
report contains a FORTRAN computer program listing. 

NEL Report 1430, - "A Numerical Analysis of HF Sector 
Ground Screen Systems," by J. M. Horn, 11 January 1967. 
Presents further numerical computations based on Wait and 
Walters 1 1963] , this time for sector screens (made up of 
radial wires) up to 128 wavelengths radius. Frequency 
10 MHz. e = 10, = .01, 1A circular impedance disc of 
12 inch square mesh of No. 10 AWG wire for various sector 
angles in 10 degree increments and observing directions 
both on the symmetry axis of a sector screen and at 10 de- 
gree increments away from this axis. It was found that if 
one stays at least 10 degrees away from the edge of a sec- 
tor then the 'gain' at any elevation angle remains within 
1 decibel of its value over a circular ^reen of the same 
radius. It was pointed out that the miu.oer of radial wires 
in a sector was doubled every time the screen radius was 
doubled. The use of such a radial wire screen can reduce 
the wire requirements 20 to 40 percent below a square mesh 
screen having the same performance. The report contains a 
FORTRAN computer program listing. 

NEL Report 1567 - "Vertical Plane Patterns of High-frequency 
Monopoles and of Elevated Vertical Dipoles with and without 
Extended Ground Planes," by J. M. Horn, (Fall 1968). Dis- 
cusses two methods of enhancing gain at low elevation angles. 
One considers the effect of long ground screens on a vertical 
A/4 monopole antenna and the other an elevated X/2 vertical 
dipole both with and without extended ground planes. Fre- 
quency 10 MHz (and a 2 to 1 frequency range above and below 
10 MHz); ground constants same as in NEL report 1359. Four 
fixed heights of the elevated dipole and ground screens of 
4 and 8 wavelengths. It is shown that an elevated vertical 
dipole at a height of A/2 with no ground plane will have the 
same low angle gain as a A/4 monopole over a 4 to 8A ground 
screen. Apparently more gain improvement is achieved by 
elevating a dipole than by constructing a ground screen, at 
least for the ground constants examined here. The improve- 
ment is greater the greater the height of the element. A 
computer program listing is available for this report. 
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The moderator invited R. G. FitzGerrell (ITS/ESSA, Boulder, 
Colo.) to comment on a method for measuring the absolute 
gain of vertically polarized HF band antennas. FitzGerrell 
then reviewed his paper [1967] where he discussed a method 
of using a horizontally polarized reference antenna as a 
gain standard, combined with an airborne platform having 
an antenna whose polarization can be switched between ver- 
tical and horizontal. Up to angles of about 30 degrees the 
'gain' of the horizontal standard is relatively insensitive 
(+ 1 dB) to the assumed ground dielectric constants and 
conductivity in the frequency range 10 to 1000 MHz. The 
results of some model measurements at 400 MHz were presented. 
J. R. Wait asked if the formula and graphs presented in- 
cluded the effect of ground on the input impedance of the 
\l2 horizontal dipole. The speaker replied in the affirma- 
tive and that his result agrees with the calculation by 
Surtees rPh.D. thesis at University of Toronto, 1.352] . 

V. Arens suggested that on the airplane shown by FitzGerrell 
the close coupling of the loops would affect the actual 
radiated power when one was rotated from a ncninal vertical 
to a nominal horizontal polarization. R. FitzGerrell re- 
plied that this effect* was corrected for experimentally in 
the data reduction of the flight test measurements. There 
was some further discussion ?mong V. Arens, R. FitzGerrell 
and H. Cottony (ITS/ESSA, Boulder, Colo.) concerning im- 
pedance variations on the aircraft antenna. 

The moderator then drew attention to the work done by A. C. 
Wilson [ 1961] at CRPL on measuring the effect of a ground 
screen of radial wires on the radiation pattern of a verti- 
cal monopole. 

The moderator suggested that Professor S. W. Maley (Univer- 
sity of Colorado, Boulder, Colo.) might like to comment on 
his work comparing measurements and the compensation theo- 
rem theory of ground screens. 

Professor Maley referred to his talk of the previous day 
in answer to this suggestion and Invited J. R. Wait to make 
any additional comments. J. R. Wait expressed the hope chat 



The actual effect referred to is the coupling of the antennas to the 
aircraft structure (Ed.). 
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the idea of a fast wave screen or screen which would exhi- 
bit negative wave tilt might (mentioned by S. Maley in 
Session IV) lead to some interesting new developments . 

W. Henry (U.S. Coast Guard, Washington, D. C.) asked how 
one can justify (in applying the compensation theorem) 
keeping the tangential H fields the same in the unperturbed 
and perturbed conditions. J. R. Wait stated that in fact 
these fieldfe aren't identical but that due to the varia- 
tional (stationary) properties of the compensation theorem 
only second order errors are made in the perturbed quantity 
being calculated. If one wishes, the compensation theorem 
can be used to verify that the perturbed cangential H 
field is only very slightly changed from its unperturbed 
value. The moderator ventured the opinion that this sta- 
tionary ptoperty of the compensation theorem was not suffi- 
ciently emphasired by its users. 

T. Kaliszewski observed that J. B. Andersen [Teleteknik, 
Copenhagen, 1965] measured the effects of an extended 
ground screen system on a log-periodic HF band antenna - 
and obtained excellent agreement with rhe predictions c.Z 
the compensation theorem. 9 " This agreement should do much 
to convince those in doubt about the credibility of the com- 
pensation theorem approach. J. R. Wait by way of reply 
suggested that it is really not necessary to rely on experi- 
ment to provide convincing evidence. In connection with a 
study of mode conversion of VLF waves at a land- sea boundary 
he showed [Wait, 1969] that a rigorous mode matching tech- 
nique and the compensation theorem approach yield the same 
first order results, provided the reflected modes are 
neglected in estimating the fields on the earth's surface 
in the later method. 

Professor R. J. King (University of Wisconsin, Madison, Wis.) 
commented that the tangential E field is oftem more sensitive 
to boundary perturbations than the H field and hence esti- 
mates of the H field are likely to be more accurate than 
those of the E field which is frequently discontinuous at 
boundaries. He then addressed some remarks to the subject 
of variable surface impedance screens, and expressed hope 
that some progress will be made on the inverse problem of 



* J. B, A derse*. Teleteknik 9, No. 2, 33-40, 1965. 
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deterraining from a specification of the desired radiation 
pattern what the required variable surface impedance must 
be. J. R. Wait then spoke briefly on some of the effects 
on the radiation pattern when a tapered impedance screen is 
used. This way it is possible to reduce the high angle 
lobing evident with abruptly terminated large screens, at 
the expense of gain at low angles. The compensation theorem 
can be used to formulate an integral equation for the im- 
pedance contrast. 

R. J. King stated that the solutions to such a problem turn 
out to be non-reciprocal and that this difficulty probably 
lies in the effect of the local ground on the impedance of 
the antennas which must be somehow properly handled in the 
formulation. (JRW s solutions referred to above were reciprocal., Ed.) 

W. L. Curtis suggested that the abrupt change of impedance 
at the edge of a screen is necessary to the attainment of 
low angle gain. 

R. D. Wengenroth was then invited to comment on some of the 
practical problems of measuring antenna system performance 
with or without ground screens. He suggested that even 
though we may some day be able to solve the synthesis prob- 
lem, it will be too difficult to realize the solution in 
practical terms because of lack of control of the environ- 
ment over a great enough distance. Another .difficulty is 
the lack of sufficiently detailed knowledge about the char- 
acteristics of probing antennas (in their environment) used 
to measure performance of ground based arrays. For HF com- 
munication systems needs the kind of results currently 
available on measured patterns and gain are usually adequate. 
In some cases however better quality data is needed and 
then the environmental influences on t^e probing antenna 
seriously limit the ability to make absolute gain measure- 
ments, but not relative gain measurements. He discussed 
the need in many cases to make measurements at very long 
distances from the nominal 'antenna' installation to take 
proper account of hills and other terrain features extend- 
ing many miles away from a rough terrain antenna site. In 
such instances a large high performance high altitude air- 
craft like the KC-135 is needed. The theoretical work to 
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determine what features of terrain are influential is a 
necessary concomitant to the experimental work. 

The moderator pointed out for the benefit of those who 
did not attend the tutorial sessions that R. H. Ott 
(ITS/ESSA, Boulder, Colo.) has developed a theoretically 
based computer executed method for predicting the pattern 
over irregular terrain, and R. Row (Sylvania, Waltham, 
Mass.) described a somewhat different approach to calcula- 
ting the field over two dimensional irregular terrain. 

V. Arens asked R. D. Wengenroth about the influence of 
the airborne loop antenna on cross polarisation measure- 
ments. R. D. Wengenroth replied that aome test results 
with both vertical and horizontal dipoles showed the 
loop to have 'good 1 (20-30 dB) rejection of the cross- 
polarized component, so that cross polarization effects 
shouldn't have corrupted the normal measurements. 

The moderator thanked the panel contributors and the 
audience in general for their thoughtful and stimulating 
discussions and called the session to a close at 3:45 p.m. 
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ADDED MATERIAL (Ed.) 

Professor S. W. Maley suggested during the round table discussions 
that laboratory models be used to study the effect of a "fast-wave" 
ground system on the radiation pattern (with reference to his paper on 
pg. 1*196), However, it was not clear how such a "fast -wave" im- 
pedance boundary is to be fabricated. 

There was also some discussion concerning the tangential field 
approximations used in applying the compensation theorem formulation 
to the se.\f -impedance calculation of a monopole. As indicated elsewhere, 
in the summarized discussions, the stationary property permits a crude 
approximation to the tangential magnetic field to yield relatively good 
self -impedance estimates. Professor R. J. King then pointed out 
formulations involving tangential magnetic fields are less influenced 
by local inhomogeneities than are representations involving the tangential 
electric fields. This point was attributed to earlier work by K. A. 
Norton (1942) on related subjects. 
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Professor R. J. King then suggested that, in the future, effort 
be directed toward the synthesis of antenna ground systems. He 
hoped to be able to synthesize the system to meet the objectives, 
e. g. , specify a variable surface impedance over the ground surface, 
and/or possibly the ground screen shape and size. 



Reference 



Norton, K. A. (1942), The polarization of downcoming ionospheric 
radio waves, FCC and NBS Rept. No. 60047, pp. 58-61 (Avail- 
able through ESSA Library, Boulder, Colorado). 
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SUMMARY OF DISCUSSIONS AND POINTS OF VIEW RAISED DURING 

REGULAR SESSIONS 



After nearly every paper, questions from the floor were raised. 
Many of these were in the nature of comments about details of the 
presented material and need not be recorded here for posterity. 
Furthermore, many of the points raised, concerning the effects of 
finite ground plane and related environmental effects, were discussed 
in the "round table" session on "Pattern Coniroi by Ground Screens. " 
The discussions and discourses at this session have been admirably 
summarized by the moderator Dr. R. V. Row. His summary appears 
on page 86 in this present volume. 

It was apparent that many authors described the results of compli- 
cated numerical calculations which arose from integral equation or 
related point -matching techniques. Since an analytical view of the 
environmental influence on antenna radiation leads to a boundary- value 
problem, it is not really surprising that various numerical techniques 
are resorted to. This trend in applied electromagnetic theory has been 
criticized by a number of prominent researchers from academia. 
Their viewpoint was elegantly put forward by Professor L. B. Felsen 
in a written communication following the conference. I am taking the 
liberty to quote from his letter as follows: 

"If I may voice again a personal concern which has to do not 
only with my session but with papers in general, there is now 
with the availability of computers the temptation to furnish de- 
tailed numerical results for special boundary value problems 
involving oddly shaped scatterers, propagation through a 
special environment, etc. Such work is justifiable if the 
numbeis are required for a special, particular application. 
However, unless such work is related to, and interpreted 
quantitatively via, approximate analytical techniques which 
grant an insight into ttie radiation and scattering mechanism, 
it does little to advance our understanding; in fact, failure 
to make such critical comparisons is to exploit inadequately 
the contribution that judiciously chosen accurate solutions 
can provide. Continued an<f. excessive emphasis on numeri- 
cal results without adequate theoretical perspective will 
tend to negate advances made through the years in propaga- 
tion and diffraction theory. I have often shuddered to think 
of the day when orientation toward numerical computation 
will be carried so far as to remove motivation for devices 
such as the Watson transformation, which not only do an in- 
efficient job much more efficiently but also add invaluable 
insight into physical phenomena. " 
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The points raised by Professor Felsen should be carefully con- 
sidered by those involved in computer-oriented approaches to 
electromagnetics. I personally share his misgivings that computer 
printout f themselves give little physical insight. On the other hand, 
iterative solutions of integral equation formulations lead to a quanti- 
tative understanding of the relevant phenomena. In particular, if the 
variational property of the compensation theorem formulation is 
utilized.an economical and accurate solution to a wide class of prob- 
lems can be obtained with ease. After all, "insight" means different 
things to different people. If all the properties of a given electro- 
magnetic environmental problem can be obtained in a compact graph- 
ical chart or nomograph, then who can say this gives less insight 
than an asymptotic result which has an uncertain region of validity 
(e, g. , it may be only valid if the range <r tends to infinity — but one 
would like to know if * is 1 meter or ) *negameter). 

Detailed Comments ^ 

Session I : Professor L. B. Felsen suggested that D. A. H^ll (author 
of paper No. 2, pg, 1-10) consider the method of Louis Cagniard in 
treating electromagnetic transient problems. (See A. T. deHocp, 
Appl. Sci. Res. Sec. B, vol. 8, 349-356, 1961.) 

Session II; Dr. J. R. Wait wondered whether R. J. Lytic? s solutions 
(in paper No. 1, pg. II- 8 ) for. the prolate spheroid could not be obtained 
by applying the reciprocity theorem to the case of a uniformly excited 
circumferential gap. The solution of this scalar problem is well known. 

Several individuals in the audience questioned the validity of Dr. 
Franceschetti's results (in paper No. 2, pg. 1-33). It was concluded 
that the results were only valid for direct current or extremely low 
frequencies. 

C. J. Sletten asked Mrs. Chatter jee (paper No. 3, pg. 1-38) how 
the boundary condition at the tip of the conical structure was satisfied. 
Dr. Wait explained that the eigenvalues were determined from che 
condition that the tangential electric fields were zero on the perfectly 
conducting spherical tip. Several questions from the floor, concerning 
the convergence of the proposed solution, were left unanswered. 



The page number in Vol. I and Vol. II where papers appear, is 
indicated. 
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Several individuals in I'.ie audience questioned Professor H. 
Kurss (paper No. 4, pg. 1-42) about the justification for his assumption 
of uniform current in the loop. Dr. Wait ventured the opinion that for 
loops near the ground, the loop circumference needed to be small com- 
pared with the skin depth in the earth. This "s much more stringent 
than merely requiring that the loop dimensions be small compared with 
a free -space wavelength. Dr. R. J. Phillips indicated that h*e and his 
colleagues at Berkeley had used a similar formulation and the paper 
was published in a recent issue of the journal Geophysical Prospecting 
(full reference not available, Ed. ). 

The interesting paper by P. R. Bannister (paper No. 5, pg. 1-45) 
provoked considerable discussion. The mathematical validity of the 
ima^ge concept was questioned out of context by some individual in 
the audience. It was stressed by the author that the results were 
valid only in the quad-static sense, i. e. , where distances were small 
compared with a free- space wavelength. 

C. J. Sletten questioned Dr. D. B. Large (paper No. 6, pg. 1-50) 
about the motivation for his finite wire calculations. The author re- 
plied that the results were used in interpreting geophysical data at 
ELF for his company. 

Session III: Several in the audience asked for clarification of R. J. 
King's analyses (papers Nos. 2 and 3, pg. 1-64, 1-66) concerning wave 
propagation over an impedance boundary. Professor King confirmed 
that his results were valid even when the refractive index of the ground 
was not large. However, the transition to a refractive index of unity 
(i. e. , free space) was not permitted. This appears to be a problem 
which has not yet been solved satisfactorily. 

Dr. Wait provided a rebuttal (in support of Arnold Sommerfeld) 
concerning the severe criticisms of the "master!' by Carson Tsao 
(paper No. 6, pg. 1-77). In the rebuttal, it was pointed out that 
Sommerf eld's results were not intended to be used in the quasi-static 
region where the height of the dipole above the ground was small com- 
pared with the wavelength. Actually, a quasi- static evaluation of 
Sommerfeld* s integrals (published in 1909!) leads to results fully 
consistent with electrostatic concepts. 

Session IV: Many of the points raised here were covered in the round 
table discussions. However, there was a question raised by Drs. Row 
and Wait about the possibility of a cross-polarized radiation from tne 
wire grid. It was concluded that the effect is small when dealing with 
a square mesh which is more or less isotropic in its reflecting 
properties. [With reference to paper No. 6, pg. 1-103, by G. A. Otteni. ] 
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Session V: In the papers (Nos. 2 and 3, pg. 1-108, 1-112) by K. K. 
Met, et al. , the rigidity boundary conditions were questioned. It was 
agreed that they are not realistic but they are convenient. 

Several people asked Dr. Oya (paper No. 5, pp. II- 39 ) how he 
managed to cope with the plasma boundary conditions on the antenna 
surface. It appears this question was avoided by assuming a current 
distribution. 

Session VI; No specific comments were recorded nor did the session 
chairman volunteer any. 

Session VII; Dr. R. V. Row questioned whether the anisotropic effects 
described by R. J. Phillips (paper No. 1, pg. 1-169) would be present 
given the very low dc magnetic fields measured near the moon. Dr. 
Phillips assured him that in the frequency range from 10-20 kHz the 
media would be anisotropic. 

E. W. Seeley (paper No. 2, pg. 1-174) was criticized for using 
efficiency to characterize the antenna gain and antenna patterns in 
Figures 5 and 6. Comparison with established antenna definitions 
was urged. 

In G. E. Webber's paper (No. 5, pg. 11-50), Dr. Wait wondered 
about the appropriateness of calling the propagation mechanism a 
surface wave or a Zenneck wave. Perhaps a guided wave or trapped 
mode might be a more accurate and descriptive term. 

Some post conference interest was Expressed in the unique phase 
dispersion method described by P. Cornille (paper No. 7, pg. 1-190). 
Actual data on a specified layered geology would be appreciated. 

Session DC; V. R. Arens asked D. C. Chang (paper No. 2, pg. 11-58) 
why so much effort was made to obtain closed form solutions rather 
than simply use a computer to evaluate the Fourier integral for all 
parameters of interest. This was answered by Dr. Chang with the 
comment that a closed form solution is preferable because of the 
greater insight into the solution which it provides and because of its 
greater utility in design. Dr. Wait (co-author) also commented later 
that the closed form solutions were needed in the numerical inversion 
of the integral equations for a finite linear antenna over a conducting 
half- space. Also, there is an economic aspect of the problem which 
is not negligible unless, of course, one has unlimited funds. 

Comments on the physical realizability of the model used by Dr. 
Curtis (paper No. 3, pg. 11-64 ) were made by Dr. Wait. He indicated 
that the postulated complex dielectric constant did not satisfy the 
Kramer s-Kronig relations and thus apparent violations of causality 
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could be expected. Dr. Wait and Dr. Curtis agreed, however, that 
the violation was not practically very important. 

A comment was made by Pitt Arnold after the session that, in 
his experimental investigations, he has never observed minima in 
the vertical radiation patterns like those shown by Dr. J. R. Wait 
(paper No. 10, pg. 11-84 ). Mr. Arnold further commented that in 
experimenting with various types of terminations for radial wire ground 
systems the radiation pattern was negligibly affected. Professor S. W. 
Maley commented that differences were probably due to different ground 
and ground system parameters. Dr. Wait later commen'ed that the 
calculated lobe structures for truncated ground systems have been 
verified by W. E. Gustafson and colleagues at the Navy Electronics 
Laboratory and W. L. Curtis and colleagues at Boeing Laboratories, 
Seattle. Of course, for a tapered system (including some radial wire 
systems), the lobe structure is less pronounced. Dr. D. B. Large 
asked how an approximate solution having no mathematical basis can 
be verified as meaningful. Dr. Wait answered that certain aspects 
of a solution such as boundary behavior or asymptotic behavior can be 
examined for correctness or approximate correctness. Also mathe- 
matical means are available for verifying that an approximate solution 
is nearly equal to the exact solution of the equation in question. 

A question concerning R. J. King's paper (No. 11, pg. 11-85 ) 
was raised about the measurement technique for the electric field 
over the dielectric wedge. Professor King explained that measure- 
ments were made by a scattering dipole technique and that great 
difficulty was experienced in synthesizing a low reflection termination 
at the edges of the dielectric wedge material. 

Mrs. G. B. Goe of ITS/ESSA made an unscheduled presentation 
on Noise and Cerenkov Radiation. (See pg. 11-114) ] 

I have prepared this summary of the discussions with the help of 
written comments supplied to me from C. J. Sletten and S. W. Maley 
for sessions which they rhaired. All the other comments, except for 
the quoted gem from Professor Felsen are based on my own (incomplete) 
recollections of the events. I accept the blame for any baises that the 
ree.de r may detect. Certainly, they do not reflect the official (or un- 
official) policy of the federal agencies of the U. S. Government which 
have sponsored this meeting. 



James R. Wait 
4 August 1969 
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CORRECTIONS TO VOLUME I 

Finite tubular antenna above a conducting half space (Paper No. 1 in 
Session HI, pp. 59-63) by D. C. Chang 

Because of errors involved in numerical computation, the scale 
factor on the abscissae in Figs. 2 and 3 should be modified by about a 
factor of 2. However, the discussion following these figures is still 
valid. The correct input conductance (in milli-mhos) of a tubular an- 
tenna of radius a = 0. 007 and h = 0. 25 X Q is tabulated in the following: 

-~- 0.27 0.3 0. 5 1.0 2.0 4.0 

X 



Sea -water 
(a =4 ohm/ m, 

e = 80 at 600 MHz) 7. 14 8.32 9.09 8. 88 8.84 8.84 

Wet-earth 

(a = 10 X 10 ~ 3 mho/m, 

e =10 at 100 MHz) 7.66 8.63 8.97 8. 86 8. 84 8.84 



When the antenna is located several wavelengths above the ground, 
the theoretical results compare favorably with experimental measure- 
ment of a single antenna in the free-space^ (8.90 milH-mhos)., For a 
more detailed discussion, the reader is referred to a paper by Chang 
and Wait (1969) listed in the references. 

t R. B. Mack, "A study of Circular Arrays, " Tech. Rept. No. 382, 
Crufts Lab. , Harvard Univ. , Cambridge, Mass., May 1, 1963. 



Reflection of waves of arbitrary polarization from a rectangular mesh 
ground screen (Paper No. 6 in Session IV, pp. 103-107) by G. A. Otteni 

On page 104, the last paragraph should read: It is assumed that. .... as 
E = Z. I, (Z. is the internal impedance of the wire and is approximately 
Z /2ff r) . . . . for an axial current I . . . .with surface impedance Z . 
In addition 

On page 107, for the sample calculations in the second paragraph, the 
grid wire spacing |a = « and b =2 m. for parallel wires. 



Numerical solution of dipole radiation in a compressible plasma with a 
vacuum sheath surrounding the antenna (Paper No. 3 in Session V, pp. 
112-116 by S.H. Lin and K. K. Mei 

The second sent nee of the last paragraph on page 114 should have 
read: This behavior of input reactance near 03 is in contrast with that 
for infinitely long cylindrical antenna found by filler (1968). 
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Phaee measurements of electromagnetic field componente (Paper No. 
7 in Session VII, pp. 190-195) by P. Coraille 

The corrected version of Fig. 4 is given below: 
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Radiation from a parallel -plate waveguide into an inhomogeneously 
filled space (Paper No. 6 in Session I, pp. 28-31} by R. J. 
Kostslnicek and R. Mittra - revised pg. 29 

for the canonical problem. The first step involves expresoing the 
fields in region B in ttrms of waveguide modes and in region £ in terms 
of a continuous spectrum (Fourier integral) representation. Similar 
transform representation are employed in regions A, D, and C. Typical 
representations for the magnetic intensity in regions B and E for a 
TM^ incident mode of unit amplitude are given by 

00 

<a y , B - cosjffi )." co S (^) e ° U> 

and 



f -£(z-£-2t) (?) 

(Hy) E = J E(a)e cos (ax) da 



o 



The longitudinal wave numbers are given by B = V(sn/b) -k in 

s o 

/' 2 2 l 12 ^ 

region B, £ = v a ~ k in E, and n - vY in C. The mode coefficients 

o o 

and spectral weight functions are B , D(a) , E(a), C(Y), and A (a). 

s 

The next step is to solve for the various mode and weight co- 
efficients by matching the transverse field components at z = it and 
z = 1 + 2t, thereby obtaining D(a) and A (a) in terms of E(a). Field 
matching is once more carried out at the z = interface for |xj £b 
and subsequently for |x| b. The resulting four equations are Fourier 
transformed and combined to yield the following relationships: 



00 

cos(Yb)TmE(Y)R'(Y) - P. V. J E'(a) f-^ - ^~-\ da, (4) 

L s 

together with two companion relationships giving the coefficients B 
and weight function C(Y) in terms of an integral relationship in- s 
volving E(a). In equations (2) and (3) we have written 
E'(a) = a sin(ab) E(a), R'(a) = R(a) exp (£Jt) and Q'(a) = Q(a) exp(-£#) 
The quotients 1/R(a) and Q(a)/R(a)are respectively the transmission 
and reflection coefficients for a uniform plane wave incident onto 
the dielectric or plasma slab at an angle given by 9 = sin" (a/k ) . 
Equations (3) and (4) represent a homogeneous integral equation ° 
for E(a) and a requirement that E(a) take on certain specified values 

e.C a = STT/b (S = 0, 1, 2, ). 



Revised pg. 30 
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The above f< — mulation is exact, and for the limiting case when 
K = 1, an exact solution for E(a) can be constructed using the 
function-theoretic technique. 



Modified Method of Solution 



The solution of (3) and (4) for E(ct) by the modified function- 
theoretic technique is accomplished by the construction of a sectionally 
holomorphic function F((j)) of a complex variable to, which has a certain 
pole-zero configuration, and specified branch singularities. 
Irtagrations in the complex plane yield results which, when compared 
with (3) and (4) and the compaion relations for B and C(y), give 
the required solutions and the normalization condition. These are: 

F(-8 ) = -b(-l) S B (l+6°)Ii (5) 
s s s s 

F(Z) = ir^(a)e jab E(a) (6) 
F(n)|i^j- + F(-n) = tttiC(Y) (7) 



ana 



F(B ) = b(-l) P B (1+6°) (8) 
P P P 



The function F{u)) may be factored into the form F(to) = F^w) T(u>) , 
where F^(ou) represents that function employed in the solution of the 
canonical problem {Mittra and Bates, 1965), and T(u)) represents the 
departure from the canonical function due to the presence of the slab. 
The preceeding integrations and comparisons together with the factored 
form of F((o) yield an auxiliary integral relationship 

T (w) - 1 + f X(z) T(z) dz . (9) 

u-8 w + z 

P o 

Equation (9) is not very convenient for numerical methods of solution 
when a is on the path a due to the singular nature of the partial 
kerntl X (z) , which has poles on o due to the surface modes excited 
with] i the dielectric slab. However, the path a may be deformed to 
say a' on which the integrand of (9) is wholly analytic. Numerical 
methods are now employed together with a process of analytic continu- 
ation, and the required values T(to) and T(-to) obtained. 

The near fields may be obtained from the modal expansion in 
region B and from the Fourier transforms in the open regions. The 
radiation fields are obtained directly from the spectral weight 
coefficients by employing tht method of saddle-point integration 
(Collin, I960). The surface modes may be obtained from the residues 
of the integrand in the transform representation in region A. 
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NOISE AND CERENKOV R> 0IATION 



G. B. Goe 

Institute for Telecommunication Sciences 
Environmental Science Services Administration 
Boulder, Colorado 80302 

Abstract 

Radiation generated by a moving electron is 
suggested as a source of noise in the ionosphere. 
The moving electron travels with a velocity which 
exceeds both the thermal velocity of the electrons 
in the ambient (homogeneous) plasma and the phase 
velocity of light in the medium in order to fulfill the 
Cerenkov condition. 

Introduction 

It has been convincingly demonstrated by using automatic gain 
control (AGC) records that the principal plasma resonances observed 
by the Alouette I topside sounder do not occur naturally but are instead 
excited by the sounder transmitter. A "noise" background is recorded, 
however, even when the transmitter is turned off. Theoretical re- 
sults suggest that the Cerenkov mechanism is useful to interpret the 
presence of the background noise and subsequent noise cut-off when 
the transmitter is turned off or the appearance of a resonance spike 
when the transmitter is turned on and the antenna is providing energy 
to the electrons. 

Theory 

Dispersion relations are derived using kinetic theory. The 
Cerenkov coherency condition is formally introduced into newly de- 
rived plasma dispersion relations in order to arrive at algebraic 
expressions for the dispersion relationships of a fully ionized 
warm, magneto-active plasma in the presence of a moving charge. 
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Comments 

The particular restriction that the moving electron be con- 
strained to travel parallel to the earth's magnetic field gives rise 
to a v/ell-defined cut-o^f for the Cerenkov zon« at the upper hybrid 
frequency (f^j I- f |j )a where fj^ is the plasma frequency and 
f y is the cyclotron frequency of the electrons. Under the same 
circumstances, the conditions near the plasma frequency are ill 
defined and, among other things, temperature dependent. 

This theoretical work was published in detail in a research 
report: Goe, G. B. , "Plasma Resonance and Radiation in the Upper 
Ionosphere, " School of Pharmacy, University of Colorado, 28 July 
1967, and has been out of print for some time. Additional copies of 
the report are, however, now available from the Department of 
Electrical Engineering, University of Colorado, Boulder, Colorado. 



